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Abstract. We prove that under certain spectral assumptions on the mon- 
odromy group, solutions of Fuchsian systems of linear equations on the Rie- 
mann sphere admit explicit global bounds on the number of their isolated 
zeros. 
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1. Introduction 

1.1. Fuchsian systems, monodromy. Consider a system of first order linear 
ordinary differential equations with rational coefficients on the complex Riemann 
sphere CP 1 . In the matrix form such a system can be written as 

X{t) = A(t)X(t), X G Mat„ x „(C), (1.1) 

where A(t) is the rational coefficients matrix and X(t) a multivalued fundamental 
matrix solution, detX(t) ^ 0. 

Assume that the system is Fuchsian, that is, the matrix-valued differential 1- 
form Q = A(t) dt has only simple poles on CP 1 , eventually including the point 
t = oo. Then the coefficients matrix A(t) has the form 

m A 

2=1 J 
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with the constant matrix residues A±, . . . , A m at the singular points t±, . . . , t m . The 
point t = oo is singular if A± + • • • + A m =/= 0, and then the residue at t = oo is the 
negative of the above sum of finite residues. 



Definition 1. The height of the rational matrix function A(t) as in (1.2), is the 
sum of norms of the residues, 

= Pi|| + --- + Pm|| + l|Aoo||, A 00 = -{A 1 + ---+A m ). (1.3) 



The singular, or polar locus £ = {t\, . . . ,t m } of the system (LI) is the union 
of all finite singular points; if t = oo is a singular point also, then we will denote 
£* = ru{oo} C CP 1 . 

Any fundamental matrix solution X{t) of the linear system ([hi]) with a rational 
matrix of coefficients is an analytic function on C ramified over the polar locus £. If 
7 is a closed loop avoiding £, then analytic continuation A 7 X(t) of the fundamental 
matrix solution X(t) along 7 produces another fundamental solution, necessarily 
of the form X(t)M 1 , where the constant invertible matrix M 7 £ Mat nx „(C) is 
called the monodromy factor. Choosing a different fundamental solution results in 
replacing M 1 by a matrix CM 7 C _1 conjugate to M 7 . 



Definition 2. The linear system (1.1) with a rational matrix A(t) and the singular 
locus £ C C, is said to satisfy the simple-loop spectral condition, if all eigenvalues 
of each monodromy factor _M 7 associated with any simple loop (non-selfintersecting 
closed Jordan curve), belong to the unit circle: 

for any simple loop 7 e 7Ti(C n £), SpecM 7 C {|A| = 1}. (1.4) 

Clearly, this condition does not depend neither on the choice of a fundamental 
solution, nor on the parametrization or even the orientation of the loop. 

1.2. Principal result. Let T C C \ £ be a simply connected domain in C \ £. 
Then any fundamental matrix solution X{t) is analytic in T and, moreover, the 
linear space spanned by all entries Xij(t), i,j = l,...,n, is independent of the 
choice of X. 

The main problem addressed in the article is to place an explicit upper bound 
on the number of isolated zeros of an arbitrary function from this linear space. 

One can rather easily see that this bound must necessarily depend on the di- 
mension n of the system and the degree m of the rational matrix function A(t). 
Simple examples suggest that the height r is also a relevant parameter. Besides, 
one should impose certain restrictions on the spectra of the residue matrices Ai 
and exclude from consideration simply connected domains that are "too spiralling" 
around one or more singularities (all these examples are discussed in details in the 



lecture notes | YakOl 



It turns out that besides those already mentioned, there are no other parameters 
that may affect the number of isolated zeros. Moreover, an upper bound can be 
explicitly computed in terms of the parameters n, m, r. Recall that an integer valued 
function of one or more integer arguments is primitive recursive, if it can be defined 
by several inductive rules, each involving induction in only one variable. This is a 



strongest form of computability, see [ Man77 , YakOl ] 



Theorem 1. There exists a primitive recursive function 9t(n, m, r) of three natural 
arguments, with the following property. 
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If the Fuchsian system (1.1) -(1.2) of dimension n x n with m finite singular 



points satisfies the simple-loop spectral condition (1.4) and its height is no greater 
than r, then any linear combination of entries of a fundamental solution in any 
triangular domain T C C \ S has no more than 9T(n, m, r) isolated roots there. 

The "existence" assertion concerning the primitive recursive counting function 
01, is constructive. The proof of Theorem [l]in fact yields an algorithm for computing 
91(n,m,r) for any input (n,m,r) S N 3 . 

It is important to stress that the bound established in Theorem |l| is uniform 
over all triangles, all possible configurations of m distinct singular points, and 
all combinations of residues Aj of total norm ^ r, provided that the simple-loop 
spectral condition holds. 

The assumption of triangularity of the domain T, as well as the focus on lin- 
ear combinations only, are not important. The general assertion, Theorem [| on 
guasialgebraicity of Picard-Vessiot fields for Fuchsian systems meeting the spectral 



condition (|LJ), will be formulated in § |2.2| after introducing all technical definitions. 



1.3. The Euler system. The simplest example of a Fuchsian system is the Euler 
system, 

X = T 1 AX, AeMat nx „(C). (1.5) 



The system (L5) can be explicitly solved: X{t) = t A = exp(Alnt). The mon- 
odromy group of the Euler system is cyclic and generated by the monodromy op- 
erator exp 2i:iA for the simple loop encircling the origin. Taking A in the Jordan 
normal form (without losing generality) allows to compute the matrix exponent 
and verify that the linear space spanned by the components of X(t), consists of 
guasipolynomials, functions of the form 

/(t) = ^t A W (lnt), A = {A 1 ,...,A„}cC, PA eCM. (1-6) 
AeA 

The spectrum A of these quasipolynomials coincides with the spectrum of the 
residue matrix A and the degree d ~ Y]a(1 + degp\) is equal to the dimension 
n of the initial system. 

Distribution of complex zeros of quasipolynomials depends on their degree and 
spectrum. For example, if A = {±i}, then the quasipolynomial f(t) = t l + t~ % — 
2 cos In t has an infinite number of positive real roots accumulating to the origin. 

On the other hand, if A C R, then at least for quasipolynomials with real coef- 
ficients, i.e., when p\ G R[lnf], the number of real positive zeros is no greater than 
d — 1 (similarly to the usual polynomials corresponding to the case A C N). 

It turns out that if A C R then not only real, but also all complex roots can be 
counted. The following result, being a particular case of Theorem |l|, will be used 
as a basis for the inductive proof of the general case as well. 



Lemma 1 (see [KY96|, Theorem 2). The number of isolated roots of any guasipoly- 
nomial ( |1.6| ) of degree ^ d with arbitrary complex coefficients but the real spectrum 
A C [— r, ?*] C R, in any triangular domain T not containing the origin t — 0, is 
explicitly bounded in terms of d and r. 



More precisely, in [KY9£] is shown that this number never exceeds Ar-\-d— 1. □ 
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1.4. Systems with pairwise distant singular points. The most difficult part of 
the proof of Theorem |l| is to treat confluent singularities, ensuring that the bounds 
on the number of zeros would remain uniform even when the distances \U — tj\ are 
arbitrarily small. If this is not the case, more precisely, if the bounds are allowed 
to depend on the configuration of the singular points, then both the formulation 
and the proof can be considerably simplified. 

The corresponding result, formulated in a proper context below (Theorem see 
7.3), differs from its unrestricted counterpart, Theorem^, by several instances. 

1. Fewer restrictions are imposed on the monodromy group: only small loops 
around singular points must satisfy the spectral condition fll.4| ). 

2. Unlike the general case, the spectral condition for such small loops can be 
immediately verified by inspection of the eigenvalues of Aj of the system: it 
is sufficient to require that all these eigenvalues must be real. 

3. Since the proof is considerably simplified, the bounds in this case are much 
less excessive, though still unlikely to be accurate. 

4. The price one has to pay is that the bounds for the number of zeros depend 
explicitly on the minimal distance between the singular points on CP 1 . These 
bounds explode as this distance tends to zero. 

This last observation on explosion of bounds notwithstanding, any three points of 
CP 1 can always be placed by an appropriate conformal isomorphism to 0, 1, oo. This 
implies computability of the number of zeros for systems having only three Fuchsian 
singular points on the sphere CP 1 , provided that the three residue matrices have 
only real eigenvalues (cf. with §1.3). Such result may be of interest for the theory 
of special functions, many of which are defined by this type of systems. 

1.5. Applications to tangential Hilbert problem. The problem on zeros of 
functions defined by Fuchsian equations, is intimately related to the tangential, or 
infinitesimal Hilbert problem [ AQ79 , Arn94 |, see | Yak01 | for the detailed exposition 
and bibliography. Recall that the problem concerns the maximal possible number 
of isolated real zeros of Abelian integrals of the form I{t) = § H=t oj, where H is 
a bivariate polynomial and lo a polynomial 1-form of a given degree d. The key 
circumstance is the fact that Abelian integrals satisfy a more or less explicitly known 
linear system of Picar d-Fuchs differential equations |AGV88| , [NY01[ 



The results formulated below, imply computability of this bound for polynomi- 
als H whose critical points are sufficiently distant from each other (Theorem ||). 
The general case still requires additional efforts to treat confluent singularities, as 
explained in [NY01|. 

Yet there is one specific hyperelliptic case when H(x, y) — y 2 + F(x), F G C[x}. 
Computability of the bound in this case under the additional assumption that all 



critical values of F are real, was proved in [NY99a| by almost the same construc- 
tion as exposed below. The difference concerns only two technical instances. The 
preliminary folding of the Fuchsian system (Appendix [X|) is made obsolete by the 
above additional assumption. On the other hand, the isomonodromic surgery ex- 
plained in and Appendix ^ in this case can be bypassed by a suitably adapted 
application of Lyashko-Looijenga theorem | Loo74 |. 



Acknowledgements. We are grateful to many people who explained us fine points 
of numerous classical results from analysis, algebra and logic that are used in the 
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2. QUASIALGEBRAICITY 

2.1. Strategy. The proof of Theorem [j] is organized as follows. First, we reformu- 
late the assertion on zeros of linear combinations of components of a fundamental 
solution, as a general claim (Theorem |^) on quasialgebraicity of Picard-Vessiot 
fields. Speaking loosely, this quasialgebraicity means algorithmic computability of 
bounds for the number of zeros of any function belonging to the function field C(X) 
generated by the entries Xij(t) of any fundamental solution X(t). 

The proof of quasialgebraicity goes essentially by induction in the number of 



finite singular points of the system (1.1). Actually, the assumptions on the system 
that are required (and reproduce themselves in the induction) slightly differ from 
those mentioned in Theorem |l|. Instead of requiring all singularities to be Fuch- 
sian, we allow one of them (at infinity) to be regular non-Fuchsian, with bounded 
negative Laurent matrix coefficients. On the contrary, we require that all finite 
singularities be aligned along the real axis and moreover belong to the segment 
[—1,1]. The section §^ contains the initial reduction, the proof of the fact that this 
additional condition can be always achieved by suitable transformations preserving 
the quasialgebraicity. 

The hard core of the inductive step is provided by the following general isomon- 
odromic reduction principle: if two given Fuchsian systems have the same mon- 
odromy in a simply connected polygonal domain U C C whose boundary is suffi- 
ciently distant from all singularities, then the Picard-Vessiot fields for these two 
systems, after restriction on U, are both quasialgebraic or both not quasialgebraic 
simultaneously. The accurate formulation of this principle is given in §||. 

The proof of the isomonodromic reduction rests upon another quite general fact, 



the bounded meandering principle [NY97, NY99t]. It provides a possibility to 



majorize explicitly the variation of argument of functions from Picard-Vessiot fields, 
along segments distant from the singular locus of a rational system. The accurate 
formulations are given in §|] while the proofs (completely independent from the rest 
of the paper) are moved to Appendix [b]. 

To carry out the inductive process, it is necessary to construct a system with 
fewer singular points, isomonodromic with a given system in a prescribed domain 
U. This isomonodromic surgery is closely related to the Hilbert 21st problem (the 



Riemann-Hilbert problem), sec [BolOO . Its analytic core is the quantitative version 
of the matrix factorization problem, see Appendix It is this step that actually 
requires extending the class of Fuchsian systems to that allowing one regular non- 
Fuchsian point (see above). 

2.2. Picard-Vessiot fields: definitions, notations. Consider a linear system 
( [Til ) with rational coefficients matrix A(t) and let X(t) be a fundamental matrix 
solution of this system. 

Consider the extension of the polynomial ring C[t] obtained by adjoining all en- 
tries Xij{t) of the chosen fundamental solution X: we will denote it by C[X] instead 
of a more accurate but certainly more cumbersome notation C[i, X\\ (t), . . . , x nn (t)]. 
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Moreover, from now on we will assume (to simplify the language) that the list of 
generators X = {Xij} of any Picard-Vessiot contains the independent variable t. 

The ring C[X] can be identified with a subring of the ring 0(t») of germs of 
analytic functions at any nonsingular point t* ^ S. The field of fractions of C[X] 
will be denoted by C(X) and can be identified with a subfield of the field A4(t*) of 
germs meromorphic at i*. The following properties of the ring C[X] and C(X) are 
obvious. 

1. The construction of C(X) does not depend on the choice of the fundamental 
solution X. For any two different nonsingular points the "realizations" of 
C(X) by fields of meromorphic germs at these points, are isomorphic to each 
other. However, this isomorphism (e.g., by analytic continuation along a path 
connecting these points, see below) is not canonical. 

2. Any element / £ C[X] (respectively, from C(X)), if identified with a germ at 

^ S, can be analytically continued as a holomorphic (resp., meromorphic) 
function along any path avoiding S. If the path is a closed loop beginning 
and ending at t», then the result of such continuation again belongs to C[X] 
(resp., C(X)). The corresponding monodromy operator is an automorphism 
of the ring (field), extending the linear transformation X i— > XM 1 on the set 
of generators of C(X). 

3. For any simply connected domain U free from singular points of the system 



(1.1), C[X] and C(X) can be identified with a subring (resp., subfield) of the 
ring 0(U) of holomorphic (resp., the field A4(U) of meromorphic) functions 
in U. These subring and subfield will be denoted respectively by C[X|J7] and 
C(X\U) and referred to as restrictions of C[X] and C(X) on the domain U. 
The field C(X) is closed by differentiation, and its subfield of constants is C. 
By definition this means that C(X) is the Picard-Vessiot extension of C(t), 



Mag94| . Somewhat surprisingly, we will not use this fact directly. 



5. The ring C[X] and the field C(X) are naturally filtered by the degree. By 
definition, degree of an element / £ C[X] is the lowest possible degree d of 
any polynomial expression Yl c k,at X a , Ck, a £ C, \a\ + k ^ d, representing 
f(t) (there may be several such combinations, as we do not assume algebraic 
independence of the entries xtj(t) of X). Degree of a fraction f = g/h, g,h £ 
C[X], f £ C(X), is the maximum of degg,deg/i for the "most economic" 
representation of /. 

6. The monodromy operators are degree-preserving. 

2.3. Computability in the Picard-Vessiot fields. To say about computability 
of bounds on the number of isolated zeros of functions from Picard-Vessiot fields, 
one has in addition to the degree of functions introduce one or several natural 



parameters characterizing the field itself (i.e., the corresponding system (1.1) for 
that matter), in terms of which the bound should be expressed. Besides, when 
counting zeros one should be aware of the multivaluedness of functions from C(X). 

The latter problem is resolved by the agreement to count zeros of functions only 
in triangular domains free from singular points. To address the former problem, we 



introduce two admissible standard classes of systems (1.1) with rational coefficients 
and in each case we list explicitly the natural parameters on which the bounds are 
allowed to depend. 

2.3.1. Fuchsian class. Let n,m,r be three natural numbers. 
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Definition 3. The Fuchsian class f(n,m,r) consists of all linear n x n-systems 
( |l.l| ) with rational matrices A(t) of the form ([L^), having no more than m singular 
points (including the one at infinity) and the height no greater than r, 

\\Ax\\ +■■■ + \\A m \\ + \\Ai + ■ ■ ■ + A m \\ < r. 

This is the basic and the most important class. Its definition is invariant by 
conformal changes of the independent variable (Mobius transformations). However, 
for technical reasons the theorem on zeros has to be formulated and proved for 
Picard-Vessiot fields built from another type of linear systems. 

2.3.2. The special class. This class consists of systems having only one regular 
eventually non-Fuchsian point at infinity. On the other hand, position of finite 
Fuchsian singularities is subject to additional restriction. 
Let again n,m,r be three natural numbers. 



Definition 4. The special class S(n, m, r) consists of linear systems (1.1 ) with the 
matrix function A(t) of the following form, 



7=1 3 fc=0 



A h A' k GMat nXn (C), (2.1) 



j 

such that: 

1. the finite singular points t±, . . . , t m are all on the real interval [—1, 1] C C, 

tj&R, Vj = l,...,m, (2.2) 

2. all coefficients Aj,A' k are real matrices, Aj,A' k £ Mat nxra (R), 

3. the total norm of Laurent coefficients of the matrix A(t) (including the sin- 
gularity at infinity), is at most r, 

m r 

£||AJ+£|K|Kr, (2.3) 

0=1 k=0 

4. the singular point t — oo is regular, and any fundamental matrix solution 
X(t) grows no faster than \t\ r as t — > oo along any ray: 

\\X{t)\\ + WX- 1 ^ < c \t\ r , c> 0, \t\ -> +oo, Argt = const . (2.4) 

The parameter r will be referred to as the height of a system from the special class. 

Remark 1. The choice of parameters may seem artificial, in particular, the as- 
sumption that r simultaneously bounds so diverse things as the residue norms, the 
growth exponent at infinity and the order of pole of A(t) at the regular singular 
point t — oo. Yet since our goal is only to prove computability without striving for 
reasonable explicit formulas for the bounds, it is more convenient to minimize the 
number of parameters, lumping together as many of them as possible. 

In both cases, the tuple of natural numbers describing the class, will be referred to 
as the parameters of the class or even simply the parameters of a particular system 
belonging to that class. Accordingly, we will refer to Picard-Vessiot fields built 
from solutions of the corresponding systems, as Picard-Vessiot fields of Fuchsian or 
special type respectively, the tuple (n, m, r) being labelled as the parameters of the 
field (dimension, number of ramification points and the height respectively). For 
the same reason we call S the ramification locus of the Picard-Vessiot field C(X). 
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2.3.3. Definition of quasialgebraicity. Let C(X) be a Picard-Vessiot field of a Fuch- 
sian or special type. Recall that by C(X\T) we denote the restriction of the field 
C(X) on a domain T C CP 1 . 

Definition 5. The field C(X) is called quasialgebraic, if the number of isolated 
zeros and poles of any function / G C(X\T) in any closed triangular domain T <g C\ 
S free from singular points of the corresponding system, is bounded by a primitive 
recursive function 9t'(<i, n, m, r) of the degree d — deg / and the parameters n, m, r 
of the field. 



Remark 2. It is important to stress that the bound for the number of zeros is 
allowed to depend only on the parameters of the class as a whole, and not on 
specific choice of the system. In other words, precisely as in Theorem |l|, in any 
assertion on quasialgebraicity of certain Picard-Vessiot fields the upper bound for 
the number of isolated zeros and poles should be uniform over: 

1. all configurations of m singular points t±, . . . ,t m , 

2. all triangles T C C \ {ii, . . . , t m }, 

3. all matrix coefficients Aj, A' k meeting the restrictions on the total norm that 
occur in the definition of the standard classes, 

4. all functions / of degree d from C(X). 

Remark 3. To avoid overstretching of the language, everywhere below outside the 
principal formulations to be "computable" is synonymous to be bounded by a prim- 
itive recursive function of the relevant parameters (usually clear from the context). 
The triangles will be always closed. 

2.4. Formulation of the results. The principal result of this paper is the follow- 
ing apparent generalization of Theorem |l|. 

Theorem 2 (quasialgebraicity of Fuchsian fields) . A Picard-Vessiot field C(X) of 



Fuchsian type is quasialgebraic if the corresponding system (1.1)— (1.2) satisfies the 
spectral condition ( |l.4| ). 

Actually this theorem is equivalent to Theorem [|. In the case m = 1 (i.e., for 
Euler systems) it follows from Lemma [l]. In the general case Theorem |^ is obtained 
as a corollary to the similar statement concerning the special class. 

Theorem 3 (main). A Picard-Vessiot field C(X) of special type is quasialgebraic 
if the corresponding system ( |1 . 1| ) , ( p.ip satisfies the spectral condition (1.4). 



2.5. Remarks. Several obvious remarks should be immediately made in connec- 
tion with these definitions and formulations. 

2.5.1. Zeros versus poles. To prove quasialgebraicity, it is sufficient to verify that 
any function / from the polynomial ring C[X] admits a computable bound for the 
number of zeros only. Indeed, for a rational function written as a ratio of two 
polynomials, the poles may occur only at the roots of the denominator, since the 
generators X are always analytic in T. 

2.5.2. Independence on generators. Quasialgebraicity of the Picard-Vessiot fields 
is essentially independent of the choice of generators. Indeed, if Y = {y±, . . . , yk} C 
C(X) is another collection of multivalued functions such that C(Y) = C(X) and 
the degrees degy^ are bounded by k G N, then any function / G C(Y) of degree d 
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in Y is a rational combination of X of degree ^ d! = 2kd. Computability in terms 
of d or dl is obviously equivalent. 

Example 1. The field C(X^ 1 ) generated by entries of the inverse matrix AT _1 (£), 
coincides with C(X). Indeed, detX € C[X] and therefore Cpf" 1 ) C C(X), and 
the role of X and X~ x is symmetric. Since A -1 is analytic outside E (has no 
poles), to verify quasialgebraicity of C(X) it is sufficient to study zeros of functions 
from C[X -1 ]. Note that in this case degA -1 = n, hence the two degrees deg x 
and deg x -i on the same field C(A) = C(A _1 ) are easily re-computable in terms 
of each other. 

2.5.3. Triangular and polygonal domains. The restriction by triangular domains 
is purely technical. As follows from Proposition [l] below, one can instead choose 
to count zeros in any semialgebraic domains in C ~ M 2 of known "complexity" 
(e.g., polygonal domains with a known number of edges). However, in this case the 
bound on zeros necessarily depends on the complexity of the domains. For instance, 
polygonal domains with a large number of edges may wind aro und sin gularities, 



spreading thus through many sheets of the Riemann surface (see [ YakOl ) . 

Proposition 1. Let S be a point set of m distinct points in M 2 and D <e R 2 \ S 
a compact simply connected semialgebraic domain bounded by finitely many real 
algebraic arcs of the degrees d\ , cfo, . . . , dk of total degree d = d± + • • • + dk- 

Then D can be subdivided into at most (d + ^)m(m + 1) + 2 connected pieces 
such that each of them lies inside a triangular domain free from points of S. 

Proof. Consider a partition of R 2 \ £ into triangles, for instance, slitting the plane 
R 2 along some of the straight line segments connecting points of S. It can be 
also considered as a spherical graph G\ with m + 1 vertices (one "at infinity" ) and 
M ^ m(m + l)/2 algebraic edges of degree 1. 

The domain D together with its complement R 2 \ D can be also considered a 
spherical graph G 2 with k vertices, k edges and 2 faces. 

Superimposing these two graphs yields a new graph with new vertices added at 
the points of intersection between the edges, and the number of edges increased 
because the newly added vertices subdivide some old edges. The number of such 
new vertices on each edge, straight or "curved" , can be immediately estimated from 
Bezout theorem. As a result, the joint graph will have at most M(d\+- ■ -+dfc+l) = 

M (d+l) straight edges, at most (Md t + 1) H h (Md k + 1) = Md + k curvilinear 

algebraic edges and at least k vertices. By the Euler formula one can place an upper 
bound on the number of faces of the joint graph, F = 2— k + E ^ 2 + M (2d +1). 

Thus D gets subdivided into a known number of connected components, each of 
them belonging to exactly one triangle T of the triangulation of R 2 \ S. This gives 
the required upper bound on the total number of triangles covering D. □ 

Having proved this proposition, we can freely pass from one simply connected 
domain to another, provided that they remain, say, bounded by a known number 
of line segments or circular arcs. 

2.6. Relative quasialgebraicity. In order to carry out the induction in the num- 
ber of singular points when proving Theorem [| we need a relative analog of quasial- 
gebraicity of a Picard-Vessiot field in a given domain U, eventually containing sin- 
gularities. This technical definition means a possibility of explicitly count zeros and 
poles not "everywhere in C" , but rather "in {/" (with the standard provision for 
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multivalued functions to be restricted on triangles). Unlike the global case, in the 
relative case the bound is allowed to depend on the distance between the boundary 
of U and the singular locus, yet this dependence must be computable. 

This distance should take into account the eventual singular point at infinity, so 
for any point set ScCwe define 

dkt(£,27Uoo) = inf m-tjl, It- 1 ]}. (2.5) 

If S is compact and disjoint with S, then this distance is strictly positive. 

Let U C C be a polygonal domain (for our purposes it would be sufficient to 
consider only rectangles and their complements). Denote by \x~\ the integer part 
of a real number x > plus 1 . 



Definition 6. A Picard-Vessiot field C(X) built from solutions of a system ( |1 . l|) 
of Fuchsian or special type, is said to be quasialgebraic in U , if: 

1. the boundary dU of the domain U is bounded and contains no singularities 



of the system fll.l|) , so that dist(<9£7, £ U oo) = p > 0, 
2. the number of isolated zeros and poles of any function / € C(A) in any 
triangle T (s U \ S free from singular points, is bounded by a primitive 
recursive function 9T"(d, n, m, r, s) of the degree d = deg/, the parameters 
n,m,r of the field and the inverse distance s = fl/p] € N. 

If the Riemann sphere CP 1 is tiled by finitely many domains U\, . . . , Uk and the 
field C(X) is quasialgebraic in each Ui, then it is globally quasialgebraic provided 
that the boundaries of the domains never pass too close to the singularities. 

3. Alignment of singular points: the first reduction 

3.1. Transformations of Picard Vessiot fields. In this section we show how 
Theorem || implies Theorems [j] and ^. 

This reduction consists in construction of an algebraic (ramified, multivalued) 
change of the independent variable z = ip(t) using some rational function ip S C(i), 
that "transforms" an arbitrary Fuchsian system (pTl|)-( |iT2| ) to another Fuchsian 
system with respect to the new variable z, having all singular points Zj on [— 1, 1] and 
only real residue matrices. Speaking loosely, one has to "fold" the Riemann sphere 
CP 1 until all singular points (both the initial singularities and the singularities 
created when folding) fall on the unit segment. Meanwhile the norms of residues of 
the obtained system should remain explicitly bounded. The general idea of using 



such transformation belongs to A. Khovanskii [Kho95|. 

The word "transforms" appears in the quotation marks since the ramified al- 
gebraic changes of the independent variable do not preserve rationality of the co- 
efficients. Indeed, if X(t) is a (germ of) fundamental matrix solution to a linear 
system (1.1) with rational matrix of coefficients AH) and z = ip(t) is a rational map 
of the t-sphere onto the z-sphere, the transform)^ X'(z) = X(tp- 1 (z)) of X(t) by 
the algebraic (non-rational) change of the independent variable ip — ip~ x does not 
in general satisfy a system of equations with rational coefficients. Performing the 
change of variables, we obtain 

dX'(z) _( ^ /dip ; 
dz 



1 The prime ' never in this article is used to denote the derivative. 
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and see that the new matrix of coefficients has in general only algebraic entries. 
This means that the held C(ip*X) generated by entries of the pullback if>*X — 
X o if; — X o tp -1 , is not a differential extension of C(z): it has to be completed. 

Remark 4. For any multivalued function / £ C(X) ramified over S, its pullback 
(transform) ip*f(z) — f(ip~ 1 (z)) by tp = tp" 1 is a multivalued function ramified 
over the union £' = tp(S) U critics consisting of the direct <p-image of S and the 
critical locus cmttp over which tp~ l is ramified. 

Lemma 2 (Alignment of singularities of Fuchsian systems). For any Picard-Vessiot 
field C(X) of Fuchsian type, there can be constructed an algebraic transformation 
ip = tp" 1 inverse to a rational map tp G C(f), and a Picard-Vessiot field C(Y), also 
of Fuchsian type, such that: 

1. i;*C(X) = C(^*X) C C(Y), 

2. for any f € C(X) of degree d, the pullback ip*f£ C(Y) has degree si d relative 
toC(Y), 

3. the ramification locus ofC(Y) consists of only finite real points, and the cor- 
responding residue matrices are all real, 

4. the degree of the map tp and the parameters of the field C(Y), in particular 
its height, are bounded by computable functions of the parameters of the field 
C(X). 

In addition, if the spectral condition ([L.4|) held for the initial Fuchsian system gen- 
erating the field C(X), then it will also hold for the Fuchsian system generating the 
field C(Y). 

3.2. Special vs. Fuchsian systems: derivation of Theorem [2| from Theo- 
rem Obviously, quasialgebraicity of C(Y), when proved, implies that its subfield 
C(ip*X) is also quasialgebraic. By Proposition [l], this means also quasialgebraicity 
of C(X). Indeed, t/3-preimages and </?-images of triangles are semialgebraic domains 
of explicitly bounded complexity, provided the degree of the map tp is known. 

As soon as all singular points of a Fuchsian system are finite and aligned along 
the real axis K, by a suitable affine transformation of the independent variable 
they can be brought to the segment [—1,1]. This transformation does not affect 
the residues of the Fuchsian system, therefore all additional conditions from the 
definition of the special class can be satisfied: in this case A' k = and the growth 
exponent at infinity is also zero. 

The folding construction shows that the question on quasialgebraicity of an ar- 
bitrary Picard-Vessiot field C(X) of Fuchsian type can be reduced using Lemma || 
to the question on quasialgebraicity of a certain auxiliary field of special type, with 
explicitly bounded parameters. 

The proof of Lemma 2 is completely independent of the rest of the paper and 
is moved to Appendix E . From this moment on we concentrate exclusively on the 
proof of Theorem ||. 

4. COMPUTABILITY OF INDEX 

4.1. Index. For a multivalued analytic function f(t) ^ ramified over a finite 
locus S, and an arbitrary piecewise-smooth oriented path 7 C C \ S avoiding this 
locus, denote by Vy(f) the index of / along 7. 

More precisely, if / has no zeros on 7, then Arg / admits selection of a continuous 
branch along 7 and the real number V 1 {f) is defined as Arg /(end) — Arg /(start) 



QUASIALGEBRAICITY OF PICARD-VESSIOT FIELDS 



13 



(this number is an integer multiple of 2ir if 7 is closed). If / has zeros on 7, then 
Vy(f) is defined as 

V~r(f) = limsupy 7 (/ - c), 

where the upper limit is taken over values c ^ f{l)- By definition, index of the 
identically zero function is set equal to along any path. 

4.2. Index of segments distant from the critical locus. Consider a linear 



system (1.1), (2.1) from the special class S(n,m,r) with the singular locus E and 
the corresponding Picard-Vessiot field C(X). 

If 7 is a bounded rectilinear segment not passing through singular points, then 



the distance p = dist(7, S U 00) from 7 to S U 00, defined as in (2.5), is positive 



Lemma 3. For any Picard-Vessiot field C(X) of special type the absolute value of 
index of any function f € C(A) along any rectilinear segment 7 C C\S is bounded 
by a computable function of d — deg /, the parameters of the field and the inverse 
distance s — [~l/p] . 

This result is a manifestation of the general bounded meandering principle, see 
|NY97 , NY99b|. It is derived in Appendix from theorems appearing in these 



eferences. 



4.3. Index of small arcs. As the segment 7 approaches the singular locus, the 
bound for the index along 7 given by Lemma [5] explodes. However, for Picard- 
Vessiot fields of Fuchsian and special types, the index along small circular arcs 
around all finite singularities remains computable. 

Lemma 4. Assume that the Picard-Vessiot field C(X) is of Fuchsian or special 
type and tj G S is a ramification point. 

Then there exist arbitrarily small circles around tj (not necessarily centered at 
tj but containing it strictly inside) with the following property. The absolute value 
of the index of any function f S C(X) of degree d along any part of any such circle 
is bounded by a computable function of d and the parameters of the field. 

In other words, index along some (certainly not all!) small circular arcs around 
every finite singularity is bounded by a computable function of the available data. 
The proof is given in Appendix |b| 

4.4. Corollaries on quasialgebraicity. The above results on computability of 
the index immediately imply some results on relative quasialgebraicity. 

Corollary 1. Any Picard-Vessiot fields of Fuchsian or special type is relatively 
quasialgebraic in any domain U free from ramification points. 

Proof. Indeed, any "polynomial" / £ C[X] has an explicitly bounded index along 
any triangle T c U by Lemma [| By the argument principle, this index bounds the 
number of zeros. The answer is given in terms of the inverse distance [dist L (dU,SU 
00)] and the parameters of the field. The assertion on zeros and poles of rational 



functions from C(X) follows from the remark in §2.5.1. □ 



Another application concerns quasialgebraicity of subfields. Assume that U is 
a polygon (say, triangle) and H = {hi, . . . , hk} C C(X\U) n M(U) is a collection 
(list) of functions from C(X) that after restriction on U are single-valued (i.e., 
meromorphic) there. Without loss of generality we may assume that hi are in fact 
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holomorphic in U\£, multiplying them if necessary by the appropriate polynomials 
from C[t] C C{X). 

Corollary 2. If C(X) is a Picard-Vessiot field of Fuchsian or special type and 
H C C(X) H 0(U) is a collection of branches holomorphic in a triangular domain 
U , then the subfield C(H) C C(X) is quasialgebraic in U . 

Proof. The same argument principle can applied to polynomials / e C[H] and the 
domain obtained by deleting from U small disks around singularities, as described 
in Lemma |J. Since / has no poles in such domain, Lemmas |^ and ^ together yield 
a computable upper bound on the number of zeros / in any triangle T C U , even 
if the sides of the triangle pass arbitrarily close to the singularities inside U. The 
case of rational functions is treated as above. □ 



Remark 5. The domain U may be not a triangle but a polygon with any apriori 
bounded number of sides. As usual, if the bound for functions / £ C(H) is to be 
expressed in terms of the degree of / with respect to the initial generators, then 
the degree of H (th e maximum of degrees of all hi) must be explicitly specified, as 



explained in §2.5.2 



5. ISOMONODROMIC REDUCTION 



The two corollaries from § 4.4 suggest that an obstruction to the (relative) quasial- 
gebraicity lies in the non-trivial monodromy of the Picard-Vessiot field. The results 
of this section show that this is the only obstruction: at least for Picard-Vessiot 
fields of special type, the monodromy group is the only factor determining whether 
the field is quasialgebraic or not. 

5.1. Isomonodromic systems, isomonodromic fields. Let U C C be a (polyg- 
onal) bounded domain and 

X = A(t)X, Y = B(t)Y, A, Be Mat„(C(f)), (5.1) 

two linear systems of the same size with rational coefficients. Denote by Sa and 
Sb their respective singular loci. 

Definition 7. The two systems are said to be isomonodromic in the domain U, if: 

1. they have the same singular points inside U, Ea H U — Sb H U = £, and no 
singular points on the boundary dU, and 

2. all monodromy operators corresponding to loops entirely belonging to U, are 
simultaneously conjugated for these two systems. 

Two Picard-Vessiot fields are said to be isomonodromic in U, if the respective linear 
systems are isomonodromic there. 

Choosing appropriate fundamental solutions (e.g., their germs at a nonsingular 
point t* € U \ £), one can assume without loss of generality that the monodromy 
factors are all equal: 

v 7 c u, x- 1 ■ A 7 V = y- 1 ■ A 7 V. 

The matrix ratio H — XY -1 of these solutions is single- valued in U (note that 
this may not be the case for another pair of fundamental solutions). Note also that 
after continuing H along a path leaving U, we arrive at a different branch of H 
which may well be ramified even in U. 
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Let U C C be a rectangle with sides parallel to the real and imaginary axes and 
having no singularities on the boundary. 

Lemma 5. Two Picard-Vessiot fields C(X) and C(Y) of the special type, isomon- 
odromic in a rectangle U , are either both quasialgebraic, or both not quasialgebraic 
in U. 

More precisely, if one of two isomonodromic fields admits a computable upper 
bound 0Ti(c?, n, m, r) for the number of zeros {poles) in U , then the other also admits 
the similar bound < yi2(d,n,m,r, s), additionally depending (also in a computable 
way) on the inverse distance s between dU and Sa U Sb U oo. 

The proof of the Lemma occupies the rest of this section. The strategy is, as- 
suming that C(X) is quasialgebraic, to prove quasialgebraicity of the compositum 
C(X,Y), the Picard-Vessiot field generated by entries of both fundamental solu- 
tions. Clearly, this is sufficient since C(Y) C C(X,Y) will in this case also be 
quasialgebraic. To count zeros of a functio n from C(X, Y ) we use a modification of 
the popular derivation-division algorithm [ IY9l| , |Rou98 ] based on the Rolle lemma 



for real functions. The role of the derivation in this modified algorithm is played 
by the operator Im of taking the imaginary part. 

5.2. Real closeness. First we observe that the Picard-Vessiot field of sp ecial type 
is closed by taking real or imaginary parts. We say that the system (|l.l|) is real, if 
the coefficients matrix A(t) of this system is real on K. 

Proposition 2. Let a d W \ S be a real segment free from singularities of a real 



system (1.1). Then for any function f £ C(X) there exist two functions from <C(X) , 
denoted respectively by Ke a f and Im^ /, such that 

Re(/| CT ) = (Re a f)\ a , Im(/| CT ) = (Im CT f)\ a . 

The degree of both functions does not exceed deg/. 

In other words, the real/imaginary part of any branch of any function / on any 
real segment a extends again on the whole set C \ S as an analytic function from 
the same field C(X) and the degree is not increased. 

Proof. Since the assertion concerns only the field, it does not depend on the choice of 
the fundamental solution X(t). Choose this solution subject to an initial condition 
X(t*) = E £ GL(n, R), where t* € a is arbitrary point on the segment. Then the 
solution remains real-valued on a, X(t) £ GL(n,R) for any t £ a, since the system 
( fLl| ) is real. The trivial operators Re CT , Im ff , defined on the generator set X of 
C(X) as 

Re CT X — X, Im CT X — 0; Re CT c = Re c, Im^ c = Im c Vc £ C, 

can be naturally and uniquely extended on the field C(X) of rational combinations 
with constant complex coefficients. They possess all the required properties. □ 

Remark 6. Note that, unlike single-valued functions, the result in general depends 
on the choice of a (this explains the notation). Yet if the restriction h £ C(X\U) 
is a single- valued branch in a domain U intersecting the real axis, then Re CT h and 
Im ff h do not depend for the choice of a segment a C U n R. This immediately 
follows from the uniqueness theorem for meromorphic functions. 
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5.3. Index and real part: the Petrov lemma. Let /: a — » C be a complex- 
valued function without zeros on a line segment a C C, and g: a — > R is the 
imaginary part of /, a real- valued function on a. We consider the case when g is 
real analytic on a and hence has only isolated zeros. The following elementary but 



quite powerful trick was conceived by G. Petrov [Pct90| 



Proposition 3. If g — hnf ^ has k isolated zeros on a, then the index of f on 
a is bounded in terms of k: 

IW)|<7r(A:+l). 

Proof. This is an "intermediate value theorem" for the map of the interval a to the 
circle, sending t G a to f(t)/\f(t)\. If V a (f) 5= tt, then the image of / should cover 
two antipodal points on the circle, and hence contain one of the two arcs connecting 
them on the sphere. But any such arc intersects both real and imaginary axis. 

The general case is obtained by subdividing a on parts along which variation of 
argument is equal to n. □ 

5.4. Compositum of Picard Vessiot fields. If C(X) and C(Y) are t wo P icard- 
Vessiot fields generated by fundamental solutions of the two systems (|5.l|) , their 
compositum C(X, Y), the field generated by entries of both fundamental solu- 
tions, is again a Picard- Vessiot field. This is obvious: the block diagonal ma- 
trix diag{X, Y} satisfies a linear system with the block diagonal coefficients matrix 
diag{A, B}. Notice that this system will be of Fuchsian (resp., special) class if both 



systems (5.1) were from this class. The dimension, number of singularities and the 
height of the system generating the compositum field, are computable in terms of 
the respective parameters of the original fields. 

As a consequence, for any function / £ C(X, Y) its index along segments or arcs 
described in Lemmas ^ and ||, is explicitly computable in terms of the degree of / 
and the parameters of the two fields. 

5.5. Demonstration of Lemma § As was already mentioned, the goal would be 
achieved if we prove that the compositum C(X, Y) in U is quasialgebraic provided 
that C(X) is quasialgebraic there. The proof begins with a series of preliminary 
simplifications. 

1. We assume that the fundamental solutions X,Y are chosen with coinciding 
monodromy factors, so that their matrix ratio H = YX^ 1 is single-valued hence 
meromorphic in U. Its degree in C(X, Y) is n + 1 (Example [l]). Thus to prove the 
Lemma, it is sufficient to prove quasialgebraicity of the field C(X, H): 

C(Y) c C(X, Y) = C{X, YX- 1 ) = C(X, H). 

Since the matrix function H has no poles in U \ S, to prove quasialgebraicity of 
C(X,H) — C(X)(H), it is sufficient to majorize the number of zeros of functions 
from C(X)[H], i.e., polynomially depending on the entries of H (cf. with Remark 



m 



2. By Remark [| Re H and ImH are unambiguously defined matrix functions, 
holomorphic in U \ £ and real on U H R. Replacing H by diag{ReiJ, ImH} if 
necessary, we may assume that H itself is real on R D U. 

3. Since H has no poles in U\S : to prove quasialgebraicity it is sufficient to con- 
sider only "polynomials in if', i.e., the ring C(X)[H] C C(X)(H) = C(X,H). By 
the assumed quasialgebraicity of C(X), the number of poles of each such function, 
represented as ^x a (t)H a , x a e C(X), is computable in terms of degXo, < degf, 
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Figure 1. Slit rectangle U 



the parameters of the field C(X) and the relative geometry of U and 2Ja U -S's- 
Therefore it remains only to majorize the number of zeros of functions from this 
ring. 

4. Any function from / = C(X)[H] can be represented as the finite sum 



/ = Y>A, n G CpO, fc< G C[ff], Im/iilR = 0. (5.2) 



=i 



The number v and the degrees degjj- Xi , deg ff /i^ are all explicitly computable in 
terms of degj^ Y f an d n - 

After all this preparatory work we can prove by induction in v, the number of 
terms in ( |5.2| ), that the number of isolated zeros of any function of such form in 
any triangle T C U \ S , is computable in terms of the parameters of the fields 
C(X) and C(Y), the degrees deg x Xi and deg H hi and the distance from U to 
Sa U U oo. The proof is very much similar to the popular "differentiation- 
division algorithm" used to majorize the number of real zeros of fewnomials in one 
variable ]IY91| , [Rou98[ , 

5. For v = 1 the assertion is immediate. Zeros of x\h\ in any triangle T C £7\ £ 
are either zeros of x\ G C(X) or zeros of ft-i. The number of zeros of the first kind 
in any triangle is bounded by virtue of quasialgebraicity of <C(X). The number 
of zeros of second kind is bounded by Corollary since hi is holomorphic with 
computable inde x alo ng the boundary. Thus the base of induction is established. 

6. For a sum (5.2) involving v > 1 terms, we first divide / by x v . This division 
may result in a loss or acquisition of a computable number of zeros or poles (oc- 
curring as poles and zeros of x v respectively) . The degree of Xj , j ^ 1, as well as 
the degree of / may increase at most by deg xi ^ d. After such division we may 
without loss of generality assume that x v = 1, that is, 

i/-i 

/ = h v + y^ Xjhj, hi, . . . , h v G C[H], X\ , . . . , Xi/ — \ 

i=i 



IS 



DMITRY NOVIKOV AND SERGEI YAKOVENKO 



Being convex, an arbitrary triangle T C U \ S may intersect R by at most one 
connected segment: TnlCcrcKxi, where a is one of the intervals between 
consecutive real singularities. Consider the domain obtained by slitting U along 
the two rays R \ a and in addition delete small disks around singular points as 
explained in Lemma [|. The result will be a simply connected domain U', still 
containing T, whose boundary consists of three types of arcs/segments (see Fig. 1): 

1. the "distant part" , the union of 6 segments belonging to the original boundary 
of dU, 

2. "small arcs" around singular points tj G £ (semi-circular or circular), 

3. at most 2m real segments <jj cl\£ (upper and lower "shores" of the slits). 
As before, the number of poles of / in U' is bounded by a computable function 

(all of these must be poles of the functions x.- L £ C(X)). Thus to majorize the 
number of zeros of / in T, it is sufficient to majorize the index of the boundary 
dU' . The contribution from the first two parts (the part sufficiently distant from 
the singular locus and small arcs around Fuchsian singularities) is computable by 
Lemmas |^ and [|. What remains is to estimate the contribution of the index of / 
along the real segments <7i with endpoints at the singular locus. 

7. The index V a {f) along each such segment a = <jj by Proposition || is bounded 
in terms of the number of zeros of Im CT /. Computation of this imaginary part yields 
(recall that x v = 1) 

v—X v—X v—X 

Im CT / = Im CT h v + } j lm a (xihi) = + K Im CT x { — 2J K^i, 

i—l i—1 i—X 

where the functions Xi = Im„ Xi again belong to C(X) by Proposition [| 

8. This computation allows to conclude that the contribution of each segment 
<7j to the total index of the boundary is majorized by the number of zeros on 
dj ■ C U \ S of an auxiliary function fj G C(X)[H] which is representable as a sum 



of ^ v — 1 terms of the form (5.2). By the inductive assumption, this number is 
bounded by a computable function of the degree, the parameters and the distance 
between U and S. 

Thus the total number of zeros of / in U' is shown to be bounded by a primi- 
tive recursive function of admissible parameters (the algorithm for computing this 



function is summarized in §5.7 below). The proof of Lemma ^| is complete. □ 



5.6. Quasialgebraicity near infinity. The singular point at t = oo is non- 
Fuchsian for systems of the special class, and hence it deserves a separate treatment. 
Let U be a complement to the rectangle, 

U = C\ {\Ret\ < 2, |Imi| < 1} 

and C(V) a Picard-Vessiot field of special type with the given parameters n,m,r 
(we retain the notation compatible with that used earlier in the proof of Lemma |s|) . 
Denote by 7 the (interior) boundary of U. By definition of the special class, all 
finite ramification points of C(Y) are outside U. 

Lem ma 6. The field C(Y) is quasialgebraic in U provided that the spectral condi- 



tion (1.4) holds for the "large loop" 7 going around infinity. 



Proof. Despite its different appearance, this is again an isomonodromic reduction 
principle: the role of the second system is played by an appropriate Euler system. 
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1. If the monodromy operator M for the large loop has all eigenvalues on the 
unit circle, then there exists an Euler system 

X = t~ 1 AX, A G Mat„(C), 

isomonodromic with the initial system in U, whose matrix residue A has only 
real eigenvalues in the interval [0,1). Since conjugacy (gauge transformation) by 
constant nondegenerate matrices does not affect the isomonodromy, one can assume 
that the norm of A is bounded, say, by n. 

2. The Picard-Vessiot field C(X) is quasialgebraic by Lemma |l| globally, hence 
in U as well. 

3. To prove quasialgebraicity of C(V) C C(X, H), where H = YX^ 1 , we follow 
the same strategy. More precisely, in this case choose U' to be the domain U slit 
along the negative real semiaxis and with a "small neighborhood of infinity" (i.e., 
complement to a disk of large radius) deleted. 

As before, we have to majorize the number of roots of functions from C(X)[7J] 
via the index of the boundary. This bound will serve all triangles T C U not 
intersecting the negative semiaxis (if T intersects the negative semiaxis, the slit 
must be along the positive semiaxis). 

4. The boundary of U' also consists of three types of segments/arcs: 

1. the "distant part" , the union of 5 segments belonging to the original boundary 
of 7 = dU, 

2. one circular "small arc" joo around infinity, 

3. 2 real segments tjj C R \ E (upper and lower "shores" of the slit). 

The only difference from the previous case is that the index of 7oo is not guaranteed 
to be computable, since Lemma ^ does not apply to the infinite singular point. Yet 
the following argument shows that for functions from C(X)[H] the index is bounded 
from above. In case the monodromy at infinity is trivial, this would correspond to 
bounding the order of pole; note that the order of multiplicity of a root (the upper 



bound for the above index) would be much harder to majorize (cf. with |Bol02|) 



Proposition 4. If the orientation on is chosen counterclockwise (as a part of 
dU'), then for any function f e C(X)[H] C C(X,Y) the index V loo (f) is bounded 
from above by a computable function of d — deg x H f and the parameter r of the 
special class only. 

Proof of the Proposition. We claim that any function / 6 C(X)[H] C C(X,Y) 
admits an asymptotic representation f(t) = i A ln fe (c + o(l)) with a real growth 
exponent A, natural k < n and a nonzero complex constant c G C near infinity 
(more precisely, in any slit neighborhood of infinity). Moreover, we claim that the 
growth exponent A is bounded from above in terms of deg x H f and the parameters 
of the field C(X,Y). 

Indeed, for functions C[X] this is obvious, since C[X] as a ring is generated over 
C by monomials t x m& t with A G [0, 1) and k ^ n. For the field C(X) it follows 
from the fact that a ratio of two functions with the specified asymptotics is again 
a function with the same asymptotical representation and the exponent being the 
difference of the growth exponents. Thus the growth exponent for / G C(X) with 
deg / = d is between — d and d. 

By definition of the special class, all entries of Y and Y~ x grow no faster than t r 
as t — > 00 along any ray. This means that H, H^ 1 may have a pole of order ^ r + 1 
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at infinity, hence the growth exponent of any function / S C[H] of degree d can be 
at most (r + l)d. 

Since the growth exponent of a sum is the maximum of the growth exponents of 
the separate terms, the growth exponent of any / S C(X)[H] is bounded in terms 
of deg / and r as asserted. 

It remains only to note that the index of any function f = t x ln fc t(c + o(l) along 
a circular arc with sufficiently large radius and the angular length less than 2ir, is 
no greater than 2ir\ + 1 . This proves the Proposition. □ 

5. The proof of Lemma || is completed exactly as the proof of Lemma Namely, 
to majorize inductively the index of any sum ( |5.2| ) along the boundary dU' , we 
estimate the contribution of the "distant part" by Lemma ^, the contribution of 
the arc 7 by Proposition ^|. The contribution of two remaining real intervals is 
majorized by the number of zeros of two auxiliary functions (the imaginary parts 
of f/xy) having the same structure (5.2) but with fewer terms, v — 1. □ 

5.7. Algorithm. The above proof of Lemma [| is algorithmic. Let / G C(Y) be an 
arbitrary function of degree d. The procedure suggested to majorize the number of 
its zeros in U (i.e., in any triangle T C U \ S) looks as follows. 

1°. Using the identity Y = H~ 1 X, we can express / as an element of C(X, H) 
and then as a ratio of two "polynomials", / = /1//2, /ij/2 £ C(X)[H]. The degree 
d! of fi in X, H is im medi ately computable, as well as the number of terms v ^ d' 
in the representation (|5.2|) . 

2°. We start from the list of two functions L = {/lj/a}, both written in the 
form of a sum ( [5.2| ). Then we transform this list, repeatedly using the following 
two rules: 

1. if / = y\ —1 Xihi is in C and x k ^ 1, then / is replaced by 2 functions, x\~ 

and f/xk = EiJi x i/ x k) + h k ; 

2. if Xk = 1, then / € C is replaced by 2m — 2 functions Im ff /, the imaginary 
parts of / on each of the real segments <7j, j = 1, . . . , 2(m — 1). 

These steps are repeated in a simple loop no more than 2d' times, after which C 
contains no more than (2m) d functions, all of them belonging to C(X). 

3°. After each cycle, the maximal degree max deg Xi of all coefficients of all 
functions in C is at most doubled (the imaginary part does not change it, division by 
the last term increases at most by the factor of 2). Thus all functions obtained after 
the loop is finished, will have explicitly bounded degrees not exceeding d" = d' ■ 2 d . 

4°. Since the field C(X) is quasialgebraic, the number of isolated zeros of all 
functions is bounded in terms of their degrees and the parameters of the field C(X). 

It remains to add the computed bounds together, adding the computable con- 
tribution coming from the index of all "distant" and "small" arcs occurring in the 
process. This contribution is bounded by a computable function of d" , the maximal 
degree ever to occur in the construction, the inverse distance from dU to the com- 
bined singular locus, and the parameters of the compositum field C(X, Y). The 
algorithm terminates. It will later be inserted as a "subroutine" into one more 
simple loop algorithm to obtain the computable bound in Theorem |3| 
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6. ISOMONODROMIC SURGERY 

Existence of a system isomonodromic to a given one in a given domain may be 
problematic in the Fuchsian class, but is always possible within the special class of 
systems. 

6.1. Matrix factorization problem. Let R — {a < \t\ < b} C C be a circular 
annulus centered at the origin on the complex plain and H (t) a matrix function, 
holomorphic and holomorphically invertible in the annulus. It is known, see [Birl3, 
GK6C| , |FM98| , |BolOO[ , that H can be factorized as follows, 

H(t) = F(t) t D G{t), D = diag{di, . . . , d n }, 

where the matrix functions F and G are holomorphic and holomorphically invertible 
in the sets U = {\t\ < b} C C and V = {\t\ > a} U {oo} C CP 1 respectively (their 
intersection is the annulus R) , and D is the diagonal matrix with the integer entries 
uniquely defined by H . 

We need a quantitative version of this result. Unfortunately, it is impossible to 
estimate the sup-norms of the matrix factors F and G and their inverses in terms 
of the matrix norm of H , as we would like to have. A weaker assertion in this spirit 
still can be proved. 

Assume that a given matrix function iJ, holomorphic and holomorphically in- 
vertible in the annulus R, admits factorization by two matrix functions F, G satisfy 
the following properties, 

H{t) = F(t)G(t), 

F, F^ 1 holomorphic and invertible in U = {|t| < b}, 

G, G _1 holomorphic and invertible in V = {a < \t\ < +00} 

and have a finite order pole at t = 00. 

A matrix function G(t), having a unique pole of finite order ^ v at infinity, can 
be represented as follows, 

G(t) = G(t) + G 1 t + G 2 t 2 + --- + G„t", Gi eMatn(C), (6.2) 

with G(t) holomorphic and bounded in {\t\ > a}. Similarly, if the inverse matrix 
G~ l (t) has a unique pole at infinity, then 

G- 1 {t) = G'(t) + G' 1 t + G' 2 t 2 + --- + G'J U , G- € Mat„(C), (6.3) 

with a holomorphic bounded term G'(t). 

We will say that the decomposition ([O]), ( |6.2D -( |6~3] ) is constrained in a smaller 
annulus R' = {a' < \t\ < b'} <e R, a < a' < b' < b, by a finite positive constant G, 
if 

WFm + wF-'m^c, \ t \<b, 

||G(t)|| + ||G / (t)||<C ) |t| > a, (6.4) 

||Gi|| + ••• + \\G V \\ + \\G[\\ + ■■■ + WG'J +u^C. 

Let R' = {a 1 < \t\ < 6'} = U' n V m R be the annulus, a' = a + |(6 - a) > a, 
V = b-Ub-a)<b, and U' = {\t\ < b'}, V ={\t\> a'}. 
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Lemma 7. Assume that the conformal width b/a — 1 and the exterior diameter of 
the annulus R are constrained by the inequalities 

b<q, b/a>l + l/q, qeN. (6.5) 

Assume that the matrix function H is bounded together with its inverse in R, 

WHm + WH^m^q', a<\t\<b, q' G N. (6.6) 

Then there exist two matrix functions F and G. holomorphic and holomorphi- 
cally invertible in smaller circular domains U' and V respectively, such that the 
decomposition ( |3.1| ) is constrained in R' by a bound C given by a computable (prim- 
itive recursive) function £(g, q') of the natural parameters q, q' . 

If H is real on intersection with the real axis R (~1 M ; then F and G can be also 
found real on R. 

The proof is given in Appendix ^[ 



6.2. Isomonodromic surgery. Consider a linear system (1.1), (2T) from the spe- 
cial class S(n, to, r) and assume that its singular locus S consists of two subsets 
S±, £2 sufficiently well apart. We prove that the system is isomonodromic to a 
"simpler" system (having fewer singular points), in appropriate simply neighbor- 
hoods of each subset Si. The proof is constructive and we show that the new 
systems can be also chosen from the special classes with the respective parameters 
(n,mi,ri), i = 1,2, explicitly bounded in terms of the parameters of the initial 
system and the size of the gap between St and Si- 



Lemma 8. Assume that the singular locus S of a system (1.1), (2.1) from the 
special class is disjoint with the annulus R = {a < \t\ < b}. 

Then there exists a system from the special class, isomonodromic with the initial 



system (1.1) in the disk {\t\ < b} and having no other finite singularities. 

The parameters of the new system are bounded by a computable function of q 
and the parameters of the initial system. 

The monodromy group of the new system satisfies the spectral condition ( |1.4[ ) if 
this condition was satisfied by the initial system ( p_.l| ). 



Proof. The construction is pretty standard, see |BolOC , AI8S]. We have only to 



verify the computability of all the relevant constraints. 

Let 7 be a loop running counterclockwise along the middle circle {\t\ = ^(a + b)} 
of the annulus and X(t) a solution defined by the condition X(t*) = E at some 
positive real point t„ inside the annulus. 

Since the matrix A(t) is explicitly bounded on 7, the Gronwall inequality implies 
that the monodromy factor M of this solution and its inverse M _1 both have the 
norm explicitly bounded in terms of the available data, hence the spectrum of M 
belongs to some annulus R* C C on the A-plane (actually we will need only the 
case when this spectrum belongs to the unit circle). 

The annulus R* can be slit along a meridian {Argi = const} so that the bound- 
ary 7* of the resulting simply connected domain will be sufficiently distant from 
the spectrum of M. Consider the matrix logarithm defined by the integral repre- 



sentation, see [Lan6£] 



A o = 7T- Q - M)' 1 InAdA e Mat„ xn (C). 
2th J„* 
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By construction, ||Ao|| is bounded by a computable function of the admissible pa- 
rameters, and the multivalued function X'(t) — t A ° has the same monodromy along 
the loop 7 as the initial system. Therefore the matrix ratio H (i) = X(t)(X' (t))^ 1 is 
single-valued, and holomorphic. By the Gronwall inequality, H, H^ 1 are bounded 
in the smaller sub-annulus R" = {a" < \t\ < b"}, a" = a+\(b-a), b" = b-j(b-a), 
together with the inverse H~ l by a computable function. 

Consider the factorization (6T) of H , described in Lemma (with R replaced by 



R"). Then the two multivalued functions Y{t) = F~ 1 {t)X{t) and Z(t) = G(t)X'(t), 
the first defined in U' and the other in V, coincide on the intersection R' = U'nV, 
so that the meromorphic matrix- valued 1-forms dY ■ Y^ 1 and dZ ■ Z -1 , defined in 
the domains {\t\ < b} and {\t\ > a}, in fact coincide on the intersection R' . Hence 
they are restrictions of some globally defined meromorphic matrix 1-form B(t)dt 
with a rational matrix function B{tY 



B(t) 



We claim that B(t) has a form (2.1) and compute (i.e., majorize) its number of 
singularities m and height r. Indeed, by construction B{t) is real on i?nR+ hence 
everywhere on R, and 

{-F-\t)-F{i)+F- l A{i)F{t), \t\<b, 
[G^G^it) + t^G^AoG-^t), \t\ > a. 

The terms F _1 (t) • F(t) and G(t) ■ G _1 (t) are holomorphic at all finite points of 
their respective domains of definition. Thus the matrix function B(t) has finite 
poles only at those poles of A(t), which are inside the disk {\t\ < a}. The residue of 
B(t) at a singular point tj is F~ 1 (tj)AjF(tj), where Aj is the residue of A(t) at this 
point. Since residues of the initial system were bounded by the parameters of the 
special class while the decomposition ( |6.l| ) was explicitly constrained, all residues 
of B(t) at all finite points are real and explicitly bounded. 

It is easy to verify that the order of pole of B(t) at t = oo and its Laurent 
coefficients are bounded: this follows from computability of the order of pole of G 



and G 1 at infinity and the explicit bounds (6.4) on the matrix Laurent coefficients 
Gi,G[. 

It remains to verify the regularity of the singular point at infinity and compute 
the growth exponent of the solution V(t). This is also obvious, since the growth 
exponents of G, G _1 are explicitly bounded whereas the growth exponents of X'(t) 
and its inverse are bounded by ||Ao||. 

The last assertion of the Lemma (on the monodromy group) is obvious: any 
simple loop avoiding the part of the spectrum inside the annulus, can be deformed 
to remain also inside the annulus, hence avoiding also singularities outside. The 
gauge transformation replacing X by does not affect the monodromy factors 

for the loops on which the transformation is defined. Hence the spectrum of all such 
loops remains the same for the initial and the constructed system. □ 

7. Demonstration of the main Theorem ||] 

In this section we prove our main result, Theorem ^ by induction in the number 
of finite singular points. The idea is to break the singular locus £ into two parts 
sufficiently apart to belong to two disjoint rectangles U\,U2, and find for each 
rectangle a system from the special class (with fewer singularities) , isomonodromic 
to the given one. Application of the isomonodromic reduction principle would yield 
then a computable global upper bound for the number of zeros of solutions. 
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7.1. Stretch and break. Given a system (1.1) from the special class ( p.l| ) with 
m finite singular points, we cover the complex plane C by three polygonal domains 
Ui,U2,Uz (two rectangles and a complement to a rectangle, all with sides parallel 
to the imaginary and real axes) so that in each of these domains the problem of 
counting zeros and poles of functions from C(X) is reduced to the same problem 
restated for the three auxiliary fields C(Yi), i — 1,2,3, each having no more than 
m — 1 ramification points. 

The first step is to stretch the independent variable so that the singular locus E 
is not collapsing to a single point. By definition of the standard class, E C [—1, 1], 
and one can always make an affine transformation (stretch) t = at' + b, a,b e R, 
< a ^ 1, \b\ ^ 1, which would send the extremal points i m i n < i m ax of E to ±1. 



Such transformation preserves the form (2.1) of the coefficients matrix, eventu- 
ally affecting only the magnitude of the Laurent coefficients A' k at infinity. However, 
because of the inequalities \a\ < 1, \b\ < 1 the impact will be limited. Indeed, the 
linear transform t = at' with \a\ < 1 can only decrease the norms ||^||, while the 
bounded translation t = t' + b, \b\ < 1, may increase them by a factor at most 2 P . 

Thus without loss of generality one can assume that the endpoints of E coincide 
with ±1. Since there are m singular points, some of the intervals (tj, ij+i) should be 
at least 2/(m — l)-long, hence the straight line I = {Ret = |(tj +tj +i )} subdivides 
the singular locus into two parts, each of them at least l/(m — l)-distant from I. 

Let U = {|Ret| ^ 2, \lmt\ < 1} be the rectangle containing inside it the 
singular locus. Denote by U\ and U2 two parts, on which U is subdivided by £, and 
let 173 be the complement to U — U\ U U2- 

For each Ei — E n Ui, i = 1,2, one can construct the annulus satisfying the 
assumptions of Lemma || with q ^ m (eventually, translated by no more than 1 
along the real axis). Thus there exist two systems, isomonodromic with the initial 



system (|l.l[ ), in the domains U\,U2, whose boundaries are at least l/(m — 1)- 
distant from all singular points. In C/3 the initial system has only one singularity (at 
infinity) and hence is isomonodromic with an appropriate Euler system generating 
the corresponding Picard-Vessiot field C(Y 3 ). 

Notice that all three auxiliary systems have fewer (at most m — 1) finite singu- 
lar points, and all belong to the appropriate special classes whose parameters are 
expressed by computable functions of the parameters of the initial system. 

7.2. Algorithm for counting zeros. Given a Picard-Vessiot field C(X) of special 
type described by the given set of parameters n,m,r and a function / e C(X) of 
degree d from this field, we proceed with the following process. 

1°. Stretch the singular locus E and subdivide C into three domains Ui,U2, Us 



as explained in §7.1, and for each subdomain construct the auxiliary field C(Y;), 
i = 1, . . . ,3, isomonodromic to C(X) in Ui. Each is again a special Picard-Vessiot 
field with no more than m — 1 points and all other parameters explicitly bounded. 

2°. By the isomonodromic reduction principle (Lemma ^ for U\, U2 and Lemma |TJ 
for C/3) the question on the number of zeros of / in Ui can be reduced to that for 
finitely many functions fi_ a from C(li), each of them having the degree (relative to 
the corresponding field) no more than 2 d . This step of the algorithm is realized by 
the "subroutine" described in § |5.7| which itself is an iterated loop. Yet the number 
of iterations is no greater than 0(2 d ). 
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3°. The steps 1° and 2° have to be repeated for each of the functions fi >a of 
known degrees in the respective Picard-Vessiot fields of special type with explicitly 
known parameters (actually, only the height has to be controlled). 

4°. After m loops of the process explained before, all functions would necessarily 
belong to the Euler fields. Their number, their degrees and the heights of the corre- 
sponding Euler fields will be bounded by a computable function of the parameters 



of the process. Since the spectral condition (1.4) for each of these fields will be 
satisfied, Lemma |l| gives an upper bound for the number of isolated zeros of each 
of these functions in any triangle free from singularities. 

5°. Assembling together all these inequalities for the number of zeros together 
with the computable contribution from the index of "distant" boundary segments 
and "small" arcs, gives an upper bound ^'(c?, n, m, r) for the number of zeros of 
the initial function / of degree d. 

6°. The running time of this algorithm is obviously bou nded by an ex plicit 



expression depending only on d,n,m,r. As is well-known | MR67 , Mac72 |, the 
bound itself must be therefore a primitive recursive function of its arguments. This 
completes the proof of Theorem H. 

7.3. Fields with pairwise distant ramification points. The arguments de- 



scribed above prove also the following result described earlier in §1.4 



Assume that the natural number s € N is so large that the singular points 



ti, . . . , t m of the Fuchsian system (El), (E2) satisfy the inequalities 



|*« — *j| < seN. (7.1) 

Theorem 4. A Picard-Vessiot field C(X) of Fuchsian type is quasialgebraic, pro- 



vided that its singular points t\, . . . ,t m satisfy the inequality (7.1) and all residues 
Ai, . . . , Am, Aoo = — Aj have only real eigenvalues. Besides the parameters 
n, m, r, the bound for the number of zeros in this case depends on the additional 
parameter s. 

Proof. To prove Theorem [|, it is sufficient to subdivide C into several rectangles 
Ui, i = 1, . . . , m each containing only one singularity, with boundaries sufficiently 
distant from the singular locus (the parameter s controls this distance from below). 
The spectral condition for the small loops jj around each singular point tj , finite or 
nor, is automatically verified: the eigenvalues n^j of the corresponding monodromy 
operators Mj G GL(n, C) are expon entials of the eigenvalues Afej of the residues 



A h fi ktj = exp2mX k ,j [BolOO , |AI88 | 



Application of the isomonodromic reduction principle allows to derive quasial- 
gebraicity of C(X) from that for each Euler field. Note that in this case there is 
no need in the isomonodromic surgery and induction in the number of singular 
points. □ 

Since any three points on the sphere CP 1 can always be placed by an appropriate 



conformal isomorphism to 0, 1, oo so that (7.1) is satisfied with s = 1, we have the 
following corollary. 

Corollary 3. Any Picard-Vessiot field with 3 ramification points, is quasialgebraic, 
provided that all three residue matrices have only real eigenvalues. 
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Appendix A. Folding of Fuchsian systems 

We first show how to complete the field C(ip*X) to a Picard-Vessiot field C(Y) 
in the case when ip(z) = \pz is the algebraic map inverse to the standard fold 
ipo : 1 1— > t 2 . Our main concern is how to place an upper bound on the height of the 
"folded" system. 

Then we construct the rational map <p as an alternating composition of simple 
folds and conformal isomorphisms of CP 1 so that it would align all singularities 



of any given Fuchsian system (1.1)— (1.2) on the real axis, while increasing in a 



controllable way the height r and the dimension n. On the final step we symmetrize 
the system to achieve the condition that all residues are real. 



A.l. Simple fold. Consider the Fuchsian system (1.1)-(JL2), assuming that both 



t = and t = oo (the critical values of the standard fold tpo) are nonsingular for it: 

m 

tj^O, j = l,...,m, J2 A J=°- C AJ ) 

i=i 

We claim that the collection of 2n functions (xi(t), . . . ,x n (t),txi(t), . . . ,tx n (t)j, 
after the change of the independent variable t = y/z satisfies a Fuchsian system of 
2n linear ordinary differential equations which can be obtained by separating the 
even and odd parts A±(-) of the matrix A{t) = A + (t 2 ) + tA-(t 2 ). 

More precisely, the (column) vector function y(z) = (x(y / z), \fz x(y / i')) satis- 
fies the system of linear ordinary differential equations with a rational matrix of 
coefficients B(z), 

dz W ' 7 ' K> 2\A+(z) A^(z) + z^E 



(A.2) 

This can be verified by direct computation. The explicit formulas (A.2) allow to 



t-A- . . A 

■ ■ i 7 



comp ute residues of B, showing that it is indeed Fuchsian under the assumptions 
(A.l). Note that even if there are two points tj — — ij with z; — Zj, the folded 
system nevertheless has only simple poles. 

The points Zj = t 2 are finite poles for B{z) with the corresponding residues Bj 
given by the formulas 

i / Aa ?- y a. 



B 3 = -( I "i 3 ), j = l,...,m. (A.3) 

In addition to these singular points (obviously, simple poles), the matrix function 
B(z) has two additional singular points z = and z = oo. The residues at these 
points are 

Inspection of these formulas immediately proves the following result. 
Proposition 5. The transform C(i/)*X) of the Picard-Vessiot field of Fuchsian 



type satisfying the condition (A.l) is a subfield of the Picard-Vessiot field C(Y) 



built from solutions of the system (A.2 
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This system has the double dimension, its singular locus on CP 1 consists of the 
squares Zj = tj, j = 1, . . . , m, and the points and oo. 



The height r' of the folded system ( A. 2 ) is bounded by a computable function of 



the height r of the initial system and the inverse distance s — [dist {0, oo})] 



[maxj = i j ... jm (|ij |, \tj |)] between the polar locus of (1.1) and the critical locus of 
the fold. 

Of course, this "computable function" can be immediately written down, but 
our policy is to avoid explicit formulas. 

A. 2. Monodromy of the folded system. We show now that the spectral con- 
dition ( |l.4D is preserved by the standard fold. This is not completely trivial. 

The proof begins by constructing a fundamental matrix solution Y(z) to the 
system flA.2| ). Let z = b be a nonsingular point for the latter, in particular, b ^ 
0,oo, and denote by t + (z), t—(z), two branches of the root *Jz for z near b, with 
t+(6) = —t-(b) denoted by a. Let X±(t), be two germs of two fundamental matrix 



solutions for the system (1.1) defined near the points ±a (they must not necessarily 
be analytic continuations of each other). Then the formal computations of the 
previous sections show that (the germ at a of) the analytic matrix function 

Y ^=Uz) + t(CU Uz)±l%)))> * ± :(C,6)-(C,±a), (A.5) 



satisfies the system (A.2). In order to show that it is nondegenerate, subtract from 
its (n + k)th row the kth row multiplied by t + {z). The result will be a block upper 
triangular matrix with two n x n-blocks X+(t+(z)) and (t-(z) — t + (z)) X-(t-(z)) 
on the diagonal, each of which is nondegenerate due to the choice of b ^ 0, oo. 

Proposition 6. If all simple monodromies of the system ( |1.1| ) have eigenvalues 
only on the unit circle, then this also holds true for the folded system (|A.2|) . 



Proof. Consider any simple loop 7 on the z-plane, starting at the point a and 
avoiding the folded locus £' = {0, z\, . . . , z m , 00}. Denote by j± its two preimages 
on the i-plane, obtained by analytic continuation of the branches t±(z) along 7. 
Two cases are to be distinguished. 

If 7 was not separating and 00, then both branches are single- valued on 7 and 
hence both j± will be closed loops on the i-plane. Denote by M± the respective 
monodromy matrices for X± along j± . Clearly, the monodromy of the whole matrix 
Y(z) along such a loop will be block-diagonal: A 1± X± = X±M±, A 7± i± = t±, 
and hence 

' X+ X_\ = (X + X_ \ (M + 
UX+ t-X- j \UX + t-X- j { M_ 



A 7 



The spectrum of such matrix is the union of spectra of M± , hence the assertion of 
the Lemma trivially holds. 

In the other case when 7 separates and 00, the branches of t± will change 
sign after going around it, hence the images j± = £±(7) will be non-closed arcs 
connecting two symmetric points ±a. However, making the second turn around 7 
will restore the closedness: both 7+ = 7- ° 7+ and "/'_ = 7+ o 7_ will be simple 
closed loops issued from the points t + (a) and t-(a) respectively. Denote by 
the monodromy matrices of the solutions X± along the corresponding loops 7^. 
Then, since A y± t± = (-l) 2 t± = t±, the square A~Y of the monodromy factor for 
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the solution Y will be block-diagonal, diag(M^, M'_) and hence its spectrum is on 
the unit circle. 

But if the eigenvalues of the square of a matrix are on the unit circle, then the 
eigenvalues of the matrix itself are also of unit modulus, and the proof of the Lemma 
is complete. □ 

A. 3. Placement of poles by a conformal isomorphism. Applying the simple 
fold that preserves the real axis can decrease the number of non-real poles, if some 
of them were on the imaginary axis. Thus combining simple folds with real shifts 
t t + c, cGM sufficiently many times, one may hope to place all poles of the 
resulting system on the real axis. 

However, it can happen that after a shift one of the poles would appear arbitrarily 
closely to t = or t = oo (or simply coincide with them). This will result in an 
arbitrarily large height of the folded system ( |A.2| ) . 

In order to correct this, instead of shifts (affine translations) between the folds 
we will make conformal isomorphisms that will keep all singularities away from the 
critical locus {0, oo} of the standard fold. Such isomorphisms can be always found. 

Let p = p m = tt/ (2m — 2) be the constant depending only on the number m of 
the points. 

Proposition 7. Let ti, . . . ,t m be any m points on CP 1 (real or not), and t\ ^ 
RP 1 . Then there exists a conformal isomorphism of the sphere CP 1 which takes 
t\ into i = +\/ — 1) preserves the real line and such that the images of the points 
t2, ■ ■ ■ ,t m are all in the annulus p m ^ \z\ ^ p" 1 . 

Proof. After an appropriate affine transformation preserving R, we may already 
assume that t\ = +i while other points tj, j = 2, . . . , m can be arbitrary. 

Using the (inverse) stereographic projection, we can identify CP 1 with the Eu- 
clidean sphere S 2 C R 3 so that the real line becomes its equator E, whereas the 
pair of points ±i is mapped into the north and south poles TV, S respectively. The 
points and oo become a pair of antipodal points on E. 

The conformal isomorphism we are looking for, will be constructed as a rigid 
rotation of the sphere around the axis NS so that some two opposite points of 
the equator, sufficiently distant from all other points t2,-->,t m , will occupy the 
positions for and oo respectively. Such a pair of opposite points always exists for 
the following simple reason. 

Consider 2m — 2 points, fe, ■ ■ ■ ,t m and their antipodes on the sphere, and a 
spherical cap of geodesic radius less than p = p m = n/(2m — 2) around each point. 
The union of these caps cannot completely cover the equator E whose length is 27T. 
Therefore there is a point b on the equator, at least p-distant from all dj and their 
antipodes. Hence both b and anti-6 are p-distant from all t%, . . . , t m . 

It remains to observe that after the rotation sending 6 to and anti-6 to oo and 
returning to the initial affine chart on CP 1 by the (direct) stereographic projection, 
the distance from to any of tj, equal to \tj\, will be no smaller than its spherical 
counterpart, and the same for the distance |l/tj| from tj to infinity. Since the 
stereographic projection of the sphere on the plane is conformal (and obviously so 
is the rigid rotation of the sphere), the result will be a conformal isomorphism as 
required. □ 

A. 4. Proof of Lemma |2|. We construct an alternate sequence of conformal iso- 
morphisms and simple folds as follows. Assuming that k ^ singularities already 
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belong to the real axis, choose one that is not in E. Using the conformal transforma- 
tion described in Proposition |t], place the chosen point to while preserving 
the real line and keeping all other singularities away from {0, 00} . Then after the 
standard fold the number of poles of the system off the real line will be by one less 
than before and the construction continues by induction in k. 

After at most m — 3 folds all singular points of the initial system will be aligned 
along the real axis (without loss of generality the first three points can be assumed 
already real). 

The chain of folds and conformal isomorphisms determines the chain of Fuchsian 
systems (defined on different copies of CP 1 ): we start from the initial system and 
then either make a conformal change of the independent variable (the procedure 
that does not change neither the residues nor the monodromy group) or replace 
the system by the folded system ( |A.2| ) of the double dimension. If the intermediate 
conformal isomorphisms are chosen as in Proposition 0, then norms of the residues 
of the suspended system after the subsequent simple folds could exceed those of the 
initial system at most by a constant factor depending only on m. All this follows 
from Proposition ||. 

The spectral condition ( |l.4| ) is also preserved in this suspension (Proposition ^J) . 
Thus the system obtained at the end satisfies all assertions of the Lemma || except 
that its residue matrices are possibly not real yet. We show that this can be 
corrected by one more doubling of dimension, this time involving reflection in the 
real axis rather than reflection in the origin. 

Together with the system x = A(t)x, whose matrix is a rational function not 
necessarily real for ( £ 1, consider its "mirror image" x* = A*(t)x* obtained by 
reflection in the real axis, denoting by A*(t) — Ait) the rational matrix function 
and x* the new dependent complex vector variable. 

Notice that A*(t) and A{t) take conjugate values for all ( e I. The variables 
u — 7}(x+x*) and v — ^(x— x*), C", together satisfy the system of equations 

'u\ = 1 / A + A* i(A-A*)\ p 
v) ~ 2 \-i(A-A*) A + A* ) \v / 

and the matrix of this system is real on R. 

Applying this last dimension doubling to the system obtained on the preceding 
step, we construct a new system with the matrix of coefficients real on M. This 
completes the proof of Lemma |[ □ 

Appendix B. Zeros of functions defined by systems of polynomial 

ORDINARY DIFFERENTIAL EQUATIONS 

In this section we derive Lemmas ^ and |] from the general the orems on zeros o f 



functions defined by systems of polynomial differential equations [ NY99b , Yak9£ ] 



B.l. A complex analog of the Vallee-Poussin theorem. Consider a linear 
nth order differential equation with analytic coefficients (t) bounded in a domain 

U c c, 

y< n > + 0l (t) y ( "- 1} + • • • + a n ^(t)y + a n (t)y = 0, 
\a k (t)\^C VteJ/cC, k = l,...,n. 
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Theorem I ( ]Yak99fl , cf. with JdlVP2S[ |). If y = /(f) ^ is a nonzero solution of 



the equation (B.l) and 7 = [toi^i] C U a rectilinear segment of length (7], then the 



index V 1 (f) of f along 7 is explicitly bounded from two sides, 
|F 7 (/)|<7r(n+l)(l + 3| 7 |C). □ 



(B.2) 



B.2. Index of small circular arcs around a Fuchsian singularity. If (B.l) 



has a singularity (pole of the coefficients) in U, the coefficients become unbounded 
and Theorem | is no t applicable. Yet if the singularity is F uchsian (i.e., if the 
coefficient ak in ( |B.l| ) has a pole of order ^ k at this point, see [AJ8S, lnc44 ), then 
the index can be majorized along small circular arcs centered at the singular point. 

Let £ — t4r be the Euler differential operator and ai(t),. . . ,a n (t) holomorphic 
germs at (C, 0). Consider the differential equation 

8 n y + ai(t)S n - 1 y + --- + a n -iSy + a n (t)y = 0, £ = t4r, t G (C, 0). (B.3) 



After expanding the powe rs £ k and division by t n the equation (B.S) becomes an 
equation of the form (BJ) with coefficients <Zfc(t) having a po le of order $C k (i.e., 
Fuchsian) at the origin t = 0. Conversely, any linear equation (B.l) having a Fuch- 
sian singular point at the origin, can be reduced to the form ( B.3 ) with appropriate 
analytic coefficients ak(t) obtained as linear combinations of holomorphic germs 
t k ak{t) with bounded constant coefficients. 

In the "logarithmic time" t — In t the Euler operator £ becomes the usual 



derivative 4z and hence the whole equation (B.3) takes the form (B.l) with the 



coefficients a^(r) = ai(e r ) defined, analytic and 27ri-periodic, hence bounded, in 
a shifted left half-plane {Rer < — B} for some B S R. Circular arcs of angular 
length ip, < ip < 2n, with center at the origin, become rectilinear segments of 
(Euclidean) length ip ^ 2-7T in the chart r . Application of Theorem | to this equation 
proves the following result. 



Corollary 4. Assume that the coefficients of the Fuchsian equation (B.3) are ex- 
plicitly bounded at the origin, |«i(0)| ^ C, i — I, . . . ,n. 

Then index V 1 (/) of any solution f ^ along any sufficiently small circular arc 
7 C {\t\ = e} of angular length less than 2-7T is explicitly bounded from two sides, 



|^(/)| <tt(» +!)(! + 6ttC). □ 



(B.4) 



A globalization of this result looks as follows. Let £ = {ti, . . . ,t m } C C be a 
finite point set, ti 5^ tj, 5 a monic polynomial with simple roots on £ and T> the 
differential operator of the following form, 

V = 8{t)±, 5(t) = (t-t 1 )---(t-t m )eC[t]. (B.5) 

Consider the linear differential equation with Fuchsian singularities on S, 



V n y + oi(t) V n ~ 1 y + ■■■ + a n _ x {t)Vy + a n (t)y = 0, 



(B.6) 



where a,k € are polynomials in t. We will refer to the sum of absolute values 
of (complex) coefficients of a polynomial p € C[t] as its height. 

Proposition 8. // the degrees a nd h eights of 6 € <C[t] and all polynomial coeffi- 
cients ak G C[t] of the equation ( B.6 ) do not exceed r e N, then for any nonzero 
solution f of this equation: 

I. if 7 is a line segment such that dist(7, S U 00) ^ 1/s for some s € N, then 
\V-y(f)\ is bounded by a computable function ofn,m,r and s; 
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2. iftj is sufficiently distant from other singularities, i.e., dist(tj, Z'xjtj}) 1/s 
for some s £ N, and ~/ £ C {\t — tj\ — e} is a simple circular arc around tj, 
then limsup e ^Q |V^ e (/)| is bounded by a computable function of n,m,r and 
s. 

Proof. Both assertions follow from Theorem | and Corollary [|. 

To prove the first assertion, one has to pass from the representation flB.6| ) to the 
expanded form 

and use the explicit lower bound for \S(t)\ on 7. To prove the second assertion 
assuming for simplicity that tj = 0, one has to substitute T> = 8(f) ■ £, where 



6 = S/t, and reduce (B.6) to the form ( JB . 3| ) . The lower bound on 5(0) coming from 
the fact that all other singular points are at least 1/s-distant from t = 0, implies 
then an upper bound for aj . □ 



B.3. Prom systems to high order equations. Consider a system (1.1) of linear 
ordinary differential equations with rational coefficients and represent its matrix 
A(t) in the following form, 

dm 

A(t) = —P(t), P = Y J Pkt\ S(t) = l[(t-tj), (B.8) 

^ ' k=a 3=1 

where P(t) is a matrix polynomial. The height of P is defined as the total norm 
of its matrix coefficients J^k II P fell- Note that the denominator S(t) is a monic 
polynomial. 



The system (1.1), (B.8) can be reduced to one high order linear equation (B.C) 



or ( B.7 ) in a variety of ways. Ye t only a few of them allow for an expli cit c ontrol of 



the height of coefficients of (BJ3) in terms of the height of P and 8 in (B.8). One of 



such reductions is outlined below (another approach, suggested in |Gri0l| , is based 
on completely different ideas and may eventually lead to better bounds). 

Let po(t,x) = c\X\ + • • • + c n x n be an arbitrary linear form on C n with con- 
stant coefficients ci, . . . ,c n £ C. Consider the sequence of polynomials p%(t, x) £ 
C[t, xi, . . . , x n ], i — 1, 2, ... , defined by the recurrent formula 



Pi+i( 



s xf ^ 9pj(t,x) ^ dpj(t,x) 
(t,x) = 5(t) — — — +2^ Qx . p ik(t)x k (B.9) 



where Pjk(t) € C[t] are entries of the polynomial matrix function P(t). In other 
words, Pi+\ is the Lie derivative of the polynomial pi £ C[i, x] along the polynomial 
vector field in C n+ , 

^)|+ E '^ illr '£- (B.10) 

j,k=i 3 

orbitally equivalent to the system (fL~l|), (B.8). Obviously, all pi are linear forms in 
the variables x\, . . . , x n . 

The chain of polynomial ideals 

(po) C (p a ,pi) C (p () ,p 1 ,p 2 ) C • • • C (p , • • ■ ,Pk) C • • • (B.ll) 
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in the Noetherian ring C[t, x] eventually stabilizes. Hence for some I < +00 one 
must necessarily have an inclusion 

Pt € (Po, ■ ■ ■ ,Pt-i), 

meaning that for an appropriate choice of polynomials b\,...,bt G C[t, x], the 
identity pi = Yli=i biPi-i must hold. Since all pi are linear forms in Xi, truncation 
of this identity to keep only linear terms implies that 

Pi(t,x) + ai(t)pe-i(t,x) H h ai(t)p (t, x) = 0, 

Oi(t) - -bi(t,0) e C[t], 

with univariate polynomial coefficients at(t) G C[t\. Restricting this identity on 
the solution x(t) of the linear system (1.1), (B.8), we immediately see that its fcth 



(B.12) 



derivative 4rkU is the restriction pk(t,x(t)) on the same solution. The conclusion 



is that the function y — p (t, x(t)) satisfies the equation (B.6) with the polynomial 
coefficients a.;(i) G C[i], having Fuchsian singularities only at the roots of S. 

The polynomial identity ( B.12[ ) is by no means unique (this concerns both the 
length £ and the choice of the coefficients Oi). However, one can construct the 
identity (B.12) of computable length t and with coefficients ai G C[i] of explicitly 
bounded heights. 

Assume that a rational matrix function A = P/S of the form ( |B.g| ) has explicitly 
bounded degree and height: 

deg P, deg 5 ^ m, height of P, 5 ^ r. (B.13) 



Theorem II (| NY99b|] , especially Appendix B). There exist explicit elementary ft 
Hons, I(n, m) and t(n, to, r) with the following properties 
For any linear system (1.1)- 



unc- 



constrained by ( B.13 ), and any linear form 
Po(x) with constant coef ficient s, one can construct polynomial identity between the 
iterated Lie derivatives (B.12) so that its length and the degrees of the coefficients 
are bounded by I and their heights by x respectively. □ 



This is a particular case of a more general theorem proved in [NY99b] for Lie 
derivatives along an arbitrary polynomial (not necessarily linear) vector field in any 
dimension. The functions [ and r can be explicitly computed and their growth rates 
for large values of n, to, r estimated. These growth rates are very large: 



\m(n,m) 



(B.14) 



as n, to 



m) ^ n m , t(n, to, r) ^ (2 - 

m(n,TO,TO.) ^ expexpexpexp(4nlnTO + 0(1)) 

The enormity of these bounds is another reason why we never attempted to 
write explicitly the bounds whose computability is asserted by our main results 
(Theorems [|, || and ||) . However, the simpler assertion of Theorem || can eventually 
be made explicit with not-too-excessive bounds based on the results of Grigoriev 
|Gri01 |. 



B.4. Proof of Lemma S. Assertion of this Lemma immediately follows from The- 
orem [l| and Proposition^ in view of the construction from the preceding section 
§ p3.3| . Indeed, a linear system from the special class S(n, m, r) has the coefficients 
matrix A(t) which, after reducing to the common denominator, can be written as 
a ratio A — P/S with explicitly bounded height and degree. 
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Any linear combination / (a function of degree 1 in the respective ring C[X] 



satisfies then the linear equation (B.6) and the first assertion of Proposition |8| 
together with the bounds provided by Theorem EL yields an upper bound for the 
index of / along any segment 7 C C \ £ sufficiently distant from the singular locus 
to admit a lower bound for \S(t)\ on it. 

To treat polynomials of higher degrees, one may cither verify by inspection that 
the proof of Theorem [n] can be modified to cover also chains starting from a poly- 
nomial po € C[t,x] of any degree d, the bounds depending on d in an explicit and 
computable way. An alternative is to notice that if x — (xi, . . . ,x n ) is solution 
of the linear system x — Ax, A — {aij(t)}fj =1 , then the pairwise products XiXj 
satisfy the linear system 



fc=l 



of dimension n(n + l)/2, whose coefficients have explicitly bounded degree and 
height. This construction can be iterated as many times as necessary to treat all 
monomials of degree ^ d. Arranging all these systems in the block diagonal matrix, 
we obtain a rational linear system Yd — Ad(t)Yd, such that the functions of degree 1 
in C[Yd] contain all polynomials of degree ^ d from C[A]. The entries of the matrix 
Ad are selected among the entries of A = A\, hence all other parameters (except 
for the dimension) will remain the same, as well and the number and location of 
the singularities. This completely reduces the polynomial case to the linear one 
treated first. 

The index of a rational function / £ C(A) of degree d does not exceed the sum 
of indices of its numerator and denominator. This observation completes the proof 
of the Lemma. □ 

B.5. Proof of Lemma The same arguments as used above (with the reference 
to the second rather than the first assertion of Proposition ||) , prove also the asser- 
tion of Lemma |] in the particular case when the system is from the special class 
and the singular point tj is finite and at least 1 / s-distant from all other finite sin- 
gularities. The bound for the index in this case depends on the additional natural 
parameter s G N, not allowed in the formulation of the Lemma. 

If we were dealing with a Fuchsian system, this would pose no problem: by a 
suitable Mobius transformation of the sphere CP 1 one could move any given point 
away from the rest of S. Since conformal changes of t preserve the Fuchsian form 
and do not affect the matrix residues, the new system will belong to the same 
Fuchsian class T{n, m, r). 

However, the special class is not invariant by Mobius transformations (and even 
affine changes of t may well affect the condition £ C [—1,1]). The arguments 
necessary involve "non-special" systems. 



Consider the system (1.1) with the coefficient matrix A(t) of the form (2.1), 
having a simple pole at the origin t\ = and make a non-affine conformal map 

z=— — ^ t = c^—, O^ceM, (B.15) 
t — c z — 1 

which preserves the origin and takes infinity to the point z = 1. 
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The new coefficient matrix A{z) = A(t)dt/dz after such transformation takes 
the form 

A m A r+2 R 

Assume that |c| < 1. Then the total norm \\B2W + • ■ ■ + ||B r +2|| of the Laurent 
coefficients Bk at the point t = 1 is explicitly bounded in terms of r (the common 



upper bound for the height of the initial system (2.1) and the order of its pole at 
infinity) . 



Indeed, the (scalar) 1-form t dt is transformed by the change of variables (B.15) 
into the rational 1-form — c k+1 (l + j^i) k = ck+luJ k, where ojk is a ratio- 

nal 1-form whose Laurent coefficients are depending only on k and bounded from 
above. This means that each matrix- valued 1-form t A' k (t) dt is transformed into a 
rational matrix 1-form with arbitrarily small matrix residues, provided that |c| > 



is sufficiently small (depending on k). In other words, the residues Bk in (B.16) 
can be made arbitrarily small by a suitable choice of c small enough 



After the transformation ( B.15| ) all finite nonzero singular points tj will occur at 



tj/(tj — c). Choosing |c| sufficiently small, one can arrange them arbitrarily 
close to the point z = 1, in particular, inside the annulus < \z\ ^ 2}. Then 
the height of the monic denominator 5(z) — (z — l) r+2 ni™( z — z j) wn ^ be then 
explicitly bounded. 

Finally, the residues Aj at all Fuchsian points remain the same, thus for all 
sufficiently small values of \c\, the transformation ( B.15| ) brings the matrix A(z) 



of the system into the form (B.16) with \\Aj\\ < r, J2k \\Bk\\ ^ 1, z i — and 
^ ^ I z j I ^ 2 for j = 2,. ..,771. After reducing it to the common denominator form 
( |B.6| ), we obtain a system with explicitly bounded height of P and (5, which has a 
singular point z = at least 1/2-distant from all other points which are still on the 
segment [—2, 2]. 

Now one can safely reduce this system to a the linear equation ( |B.6| ) as described 
in Theorem || and apply the second assertion of Proposition ^ to obtain computable 
two-sided bounds on the index of small circular arcs on the z-plane around the 
origin. The images of these arcs on the t-plane will be arbitrarily small circular 
arcs around tj, though not necessarily centered at this point. The proof of Lemma ^ 
is complete. □ 

Appendix C. Quantitative factorization of matrix functions 

IN AN ANNULUS 

As often happened before, the proof with quantitative bounds is essentially ob- 
tained by inspection of existing qualitative proofs, supplying quantitative argu- 
ments when necessary. 

Th e stand ard proof of matrix factorization theorem that we "quantize" , is that 
from |BolOC |. Namely, for the matrix function H(t) sufficiently close to the (con- 



stant) unity matrix E, one can prove existence of holomorphic and holomorphically 
invertible decomposition H = FG, including the point at infinity (i.e., with both 
G and G~ 1 bounded in the outer domain). Then, using approximation by rational 
functions, we prove that there exist meromorphic and meromorphically invertible 
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decomposition with explicit bounds on the number of zeros and poles of the de- 
terminants det F, det G. Finally, we show that all zeros and poles can be forced to 
"migrate" to the single point at infinity, while retaining control over all magnitudes. 

Throughout this section R is the annulus {a < \t\ < b} with a, b constrained by 
the natural parameter q as in ( |6.5[ ), and U, V denote respectively the interior and 
exterior disks, U = {\t\ < b}, V = {\t\ > a} C C, so that R = U n V. 

C.l. Holomorphic factorization of near-identity matrix functions. 

Lemma 9. There exists a computable function 91(g) of one integer argument q 
with the following property. 

For a give n an nulus R whose dimensions are determined by the integer parameter 
q € N as in ( |6.5| ), and any holomorphic invertible matrix function H in the annulus, 
close to identity enough to satisfy the condition 

\\H(t) - E\\ < 1/Vl(q), teR, (C.l) 

there exists the decomposition H = FG with F and G holomorphic and holomorphi- 
cally invertible in the respective circular domains U and V . The matrix functions 
F, G satisfy the constraints 

\\F(t)\\ + \\F-\t)\\ < K(q), \\G(t)\\ + WG-'m < m(q) (C.2) 

in their domains. 

Scheme of the proof. The identity H = FG can be considered as a functional equat- 
ion to be solved with respect to F and G for the given matrix function H . Repre- 
senting the three functions as H = E + C, F = E + A, G = E + B, this identity 
can be rewritten as 

C = A + B + AB. (C.3) 

This equation is nonlinear with respect to the pair of matrix functions (A, B). Its 
"linearization" is obtained by keeping only "first order" terms, 

C = A + B. (C.4) 

Solution of this linearized equation can be immediately given by the Cauchy inte- 
gral. The boundary 7 of the annulus R consists of two circular arcs, the interior 
arc 7 a and the exterior arc 7^ (properly oriented). Writing 

1 £C{z)dz _ 1 I C(z)dz ( I C(z)dz 



2ni 71 z — t 2iri L z — t z — t 

" / "la " \b 

The first term which we denote by B(t), is holomorphic in V = {\t\ > a}. The 
second term is respectively holomorphic in U = {\t\ < 6} and after being denoted by 
A(t) , constitutes together with B a holomorphic solution of the linearized problem 



( |C.4D . This solution considered as an integral operator C(-) 1— * (A(-),B(-)) is of 
explicitly bounded norm. Indeed, if || C{f) \\ is bounded on 7 by 1, then for each t e R 
at least one of the distances dist(£,7 a ) or dist(£,7b) will be no smaller than l/2p, 
where q is the width parameter of the annulus as described in ( |6.5| ). This means that 
the norm of one of the integrals is explicitly bounded in terms of <?, but the second 
integral is then also bounded by virtue of ( p.4[ ). Note that this boundedness is 
specific for the annulus having two disjoint parts 7b of the boundary 7 (normally 
the Cauchy operators have unbounded norms). 



The nonlinear equation (C.3) can be now solved by the Newton method. If 



H = E + C and ||C|| ^ e, then factorization of H can be reduced to factorization 
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of (E + Ay 1 H(E + B)^ 1 : if the latter is factorized as F'G' as required, the 
former can also be factorized as (E + A)F' ■ G'(E + B). The new difference C" = 
(E + A)-~>-(E + C){E + B)- 1 - E is of quadratic order of magnitude, |jC"|| Le 2 , 
where L — L(q) is the norm of the above described integral operator solving the 
linearized equation. The process of linearization and solving the linearized equation 
can be continued. 

The super-exponential convergence of Newton-type iterations guarantees that if 
the original difference max f6 j{ H-ff(i) — E\\ wa s smaller than, say, l/2L(q), then the 
process converges and the ultimate result, constructed as the infinite product of 
matrices fast converging to the identity E in the respective domains U, V, will be 
of matrix norms explicitly bounded in terms of q only. The accurate estimates can 



be extracted from [BolOO if necessary. 



If C is real on R n R, then clearly the solution of the linearized equation (C.4) 
is also real on the real line, and this condition is preserved in the iterations of the 
Newton method. This completes the proof of the Lemma. □ 

C.2. Meromorphic solutions of the factorization problem. The following 
construction is also standard. 

Proposition 9. Any matrix function H satisfying the assumptions of Lemma [J 
can be approximated by a rational matrix function 

d 

M(t) = M kt\ M k £ Mat„(C), (C.5) 
k=-d 
r-i/ 



so that the product M(t)H (t) satisfies the condition (C.l): 

\\M(t)H-\t) - E\\ < l/m(q) 

in the smaller annulus R' = {a' < \t\ < b'} (e R, a' = a+ \ {b— a) . b' = b— j(b — a). 

The degree d and norms of the matrix coefficients Mh will be bounded by primitive 
recursive functions of q and q' . If H is real on R D M, then M will also be real on 
R. 

Proof. Consider the Laurent expansion for H , converging uniformly on R' d R since 
H is holomorphic in the annulus R. The rate of this convergence can be explicitly 



estimated in terms of q' and q, see (3.5) and (|6.6| ). Thus one can immediately 



estimate the number of terms of the Laurent expansion that would be sufficient to 
keep to satisfy the required accuracy of the approximation so that 

\\H(t)-M(t)\\ < l/q'm(q). 

Then H~ 1 AI — E and MH^ 1 — E will be as small as asserted on R' . The rest is 
obvious. □ 

This Proposition together with Lemma ^ applied to i/ _1 M, guarantees for an 
arbitrary holomorphically invertible matrix function H existence of a factorization 
of the form 

H = FG, G = G'M, 

, , (C6) 
F, G holomorphic in U , V respectively, M rational, 

where F, G' are holomorphic, holomorphically bounded and invertible in U' = {\t\ < 
b'} g U, V = {\t\ > a'} <s V, and are bounded together with their inverses there 
by a computable function of q, q' only. 
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C.3. Expulsion of poles of G 1 . Proof of Lemma [7]. The term F in the 



decomposition (JCJj) already meets the requirements imposed in (3.4). 

To achieve the conditions required from the second term G, we multiply it from 
the left by several rational matrix functions, at the same time multiplying F from 



the right by their inverses. The construction, due to G. Birkhoff [Birl3|, is utterly 
classical. 

The matrix function G has a unique pole of order ^ d at infinity, but may have 
a number of " zeros''' (poles of the inverse matrix G _1 ). Their number and position 
are determined by zeros of det G — det G' det M, that is, by zeros of the rational 
function det M(t) of known degree, since the term G' is holomorphically invertible. 
Thus G _1 has no more than a computable number of poles z\,...,Zk, k being 
bounded by a computable function of q, q' . We are interested only in those poles 
that are in V'. 

Let t = z be any root of the determinant detM(t). The matrix G(z) has non- 
trivial left null space, since det G(z) — 0. Consider the constant matrix C = C z 
built as follows: its first row is the (left) null vector of G(z) of unit Hermitian norm, 
while the other rows are chosen among the rows of the identity matrix E in such 
a way that both ||C|| and ||C _1 || arc bounded by n. This is always possible: it 
is sufficient to delete from the coordinate vectors in C™ the vector which has the 
largest Hermitian scalar product with the null vector. 

By construction, the product C z G(t) at the point t = z has the zero first row. 
Hence the product diag{r z (i), 1, . . . , 1} • CG(t), r z (t) = \z\/ (t — z), continues to be 
holomorphic at the point z. The correction term Q z (t) = diag{r z (i), 1, . . . , 1} • C is 
bounded together with its inverse in the thinner annulus R" = {a" < \t\ < b"} <s R', 
the "middle belt" of R', a" = a!+\{b'-a% b" = b' -\{b' -a'). Note that this bound 
is uniform as z varies in V = {b' < \z\}. Besides, Q z (t) is bounded together with 
its inverse also in U'. All these bounds are computable functions of q, q'. Replacing 
the decomposition H = FG by H — FQ~ r Q z G, we obtain a new factorization of 
H that has fewer poles of the right term in V (both terms remain holomorphic and 
the left is holomorphically invertible in U'). 

This expulsion of a pole of G _1 from t = z to t = 00 has to be repeated for each 
pole Zj <G V' of the initial matrix G = G'M (the order is inessential). As a result, 
we replace the decomposition ( p.6[ ) by 

H = FQ^ 1 ■ QG, Q = Q(t) rational without poles or zeros in R' . (C.7) 

By construction, both QG and (QG)^ 1 have only one (multiple) pole at t — 00, 
and all four matrices FQ~ l , QG and their inverses, are explicitly bounded in R" 
in the sense of the norm by a computable function of q, q' . 

The order of pole of QG at infinity is no greater than the initial order d of pole 
of M. The order of pole of (QG)^ 1 does not exceed k, the degree of detM (the 
number of singular points that were expelled) . Finally, the norms of the Laurent 
coefficients of (QG)^ 1 at infinity can be majorized using the Cauchy estimates 
(differentiating the Cauchy integral along the inner boundary of the annulus R"). 
These bounds will prove the assertion of Lemma 0. 

The only remaining problem is to choose Q being real on K provided that G is 
real there. In this is the case, then non-real "zeros" of G come in conjugate pairs 
z, z and the corresponding left null vectors (rows) C™ are complex conjugate. 

In this case both zeros can be expelled if the corresponding rational matrix factor 
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Qz,z{t) has the first row of the form 

M 



t - 



t e 



This vector function has the norm bounded from above away from the pair of points 
{z, z} uniformly over all z £ V, and as before the inverse Q^ 1 is of bounded norm 
provided that the other rows of Q were appropriately chosen among the rows of the 
identity matrix E exactly as before. 

This remark concludes the proof of Lemma |7[ □ 
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